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ABSTRACT 


In  the  design  and  analysis  of  a  digital  satellite 
communications  system  there  is  an  acute  need  to  have  the 
ability  to  predict  the  system's  bit  error  rate  (BER) . 

After  defining  the  ideal  performance  prediction  tool,  those 
that  have  been  proposed  or  are  now  in  use  are  examined  and 
all  are  seen  to  fit  into  one  of  two  general  categories: 

The  analytical  methods  and  the  Monte  Carlo  simulation- 
based  approaches.  A  comparative  review  of  these  predictors 
is  presented.  All  are  found  to  fall  short  of  the  ideal 
although  the  simulation-based  estimators  are  much  closer  than 
any  of  the  analytical  ones.  The  best  performance  prediction 
tools  are  seen  to  be  an  application  of  importance  sampling 
to  Monte  Carlo  simulation  proposed  by  Shanmugam  and  Balaban 
and  a  graphical  extrapolation  method  proposed  by  Weinstein 
that  entails  linearization  of  the  tail  of  the  output 
signal's  pdf.  A  new  predictor  is  proposed  that  combines 
the  techniques  of  Shanmugam  and  Balaban  with  those  of 
Weinstein  in  an  attempt  to  utilize  the  advantages  of  each. 
This  Modified  Extrapolation  technique  is  then  evaluated  and 
found  to  be  better  than  Weinstein's  method  but  just  slightly 
worse  than  the  importance  sampling  approach  for  a  system 
BER  =  10  6  with  a  normalized  error  =  0.1.  At  smaller  BERs, 
the  structure  of  the  proposed  technique  indicates  that  it 
will  again  be  superior  to  Weinstein's  approach  and  that  it 
will  even  become  better  than  Shanmugam  and  Balaban 's  method. 
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1.0 


INTRODUCTION 


A  typical  digital  satellite  communications  system 
consists  of  two  earth  stations  and  a  channel  containing  an 
uplink,  a  downlink,  and  the  satellite  transponder.  In  the 
process  of  designing  such  a  system,  it  is  absolutely  essential 
that  the  designer  have  the  capability  to  predict  its 
performance  under  conditions  similar  to  those  that  will  be 
encountered  during  actual  operation.  Having  a  performance 
prediction  tool  allows  the  comparisons  of  various  filters, 
modulation  schemes,  amplifier  operating  points,  etc.,  in 
the  quest  for  a  "best"  system  to  fit  the  design  criteria. 
Further,  the  tool  permits  the  determination  of  a  system’s 
sensitivity  to  perturbations  such  as  phase  imbalances,  linear 
and  nonlinear  (AM-AM,  AM-PM)  distortions,  and  time,  phase, 
and  frequency  jitter. 

While  predicting  the  performance  of  the  typical  satellite 
system  involves  the  determination  of  several  different 
performance  indices,  perhaps  none  is  quite  so  useful  and 
important  as  the  symbol,  or  bit,  error  probability  [P(e)]. 
Commonly  known  as  the  bit  error  rate  (BER) ,  this  parameter 
is  the  prevalent  one  used  to  make  qualitative  and  quantitative 
evaluations  of  the  performance  of  digital  communications 
systems.  When  operation  is  over  a  linear,  Gaussian  channel 
with  no  bandwidth  restrictions,  calculation  of  the  error 
probability  of  a  system  is  a  very  straightforward  and  easy 
task  [1] . 

Performance  prediction,  however,  is  an  attempt  to 
determine  how  well  a  realistic  system  will  perform.  In 
digital  satellite  communications  this  means  taking  into 
consideration  the  facts  that  the  channels  used  are  frequently 


bandlixnited  and  often  nonlinear,  as  is  the  case  when  a 
travelling-wave  tube  amplifier  (TWTA)  is  operated  near  its 
saturation  point.  The  bandlimiting  causes  the  modulated 
pulse  to  suffer  some  spreading  across  time  thus  introducing 
intersymbol  interference  (ISI) .  In  addition,  the  channel 
noise  is  occasionally  non-Gaussian.  With  a  system  having 
any  of  these  characteristics,  the  task  of  calculating  the 
BER  performance  is  quite  difficult. 

Finding  a  solution  to  this  problem  is  a  matter  of 
finding,  or  creating,  an  acceptable  prediction  tool.  In 
order  to  meet  the  requirements  stated  earlier,  the  tool  must 
be  very  flexible  with  regards  to  system  configurations  and 
components.  On  the  other  hand,  it  is  obvious  that  the  tool 
should  produce  accurate  results  without  an  unreasonable 
computational  and  analytical  burden.  Unfortunately,  these 
two  requirements  have  proven  to  be  unreconcilable  [2].  Out 
of  these  requirements  have  grown  the  two  techniques  that 
form  the  basis  for  all  performance  prediction  tools.  Monte 
Carlo  simulation  is  an  approach  in  which  the  system  is  modeled 
in  a  component  by  component  (building  blocks)  manner  with  the 
simulation  being  run  in  the  time  domain  and  the  noise  treated 
as  a  vector  of  random  samples.  This  can  be  used  to  model 
the  performance  of  any  system,  with  any  configuration,  to 
any  desired  accuracy.  However,  due  to  the  extremely  low  bit 
error  rates  seen  in  satellite  systems,  typically  on  the  order 
of  10“6,  the  time  involved  to  produce  reliable  results  is 
often  unacceptable  because  of  large  sample  size  requirements. 
On  the  other  hand  there  are  purely  analytical  methods  that 
can  quickly  produce  acceptably  accurate  results  using  one  or 
more  integrated  mathematical  models  of  the  system.  However, 
the  assumptions  needed  to  derive  analytically  tractable 
results  are  often  too  restrictive,  thus  limiting  application 
to  a  few  easily  characterized  systems. 
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It  is  well  known  that  the  primary  causes  of  detection 
error  in  digital  communication  are  thermal  noise,  typically 
additive  white  Gaussian  noise  (AWGN) ,  and  the  intersymbol 
interference  that  arises  in  the  communications  channel. 
Calculating  the  error  probability  for  a  satellite  system, 
therefore,  is  equivalent  to  finding  a  prediction  tool  that 
will  account  for  the  effects  of  noise  and  ISI  on  the  signal. 

This  means  taking  into  consideration  the  effects  of  the 
nonlinearly  transformed  uplink  noise  and  ISI  in  addition 
to  those  of  the  ISI  and  noise  arising  on  the  downlink. 

Further,  the  tool  must  all  the  while  incorporate  the  conflicting 
requirements  of  speed  and  flexibility  mentioned  previously. 

Over  the  past  two  decades,  a  considerable  amount  of  work 
has  been  done  towards  solving  the  various  parts  of  this 
problem.  There  have  been  analytical  approaches  [3]-[30], [47] , 
approaches  based  on  Monte  Carlo  simulations  [ 31] — [ 46 ] ,  and 
even  a  hybrid  approach,  but  as  of  yet,  there  does  not  seem 
to  be  a  clear  understanding  of  the  overall  effect  that  uplink 
noise  and  intersymbol  interference  preceding  a  nonlinearity 
have  on  the  BER. 

Almost  all  of  the  early  work,  and  the  vast  majority  of 
what  has  been  done  so  far,  can  be  placed  into  the  realm  of 
the  analytical  techniques.  A  great  amount  of  effort  has  been 
given  to  the  study  of  the  effects  of  ISI  on  the  BER  performance 
of  a  system  but,  in  general,  either  a  linear  channel  has 
been  used  [3]- [13]  or  if  a  nonlinear  channel  has  been  assumed, 
then  the  uplink  noise  has  been  either  ignored  or  considered 
to  somehow  bypass  the  nonlinear  transformation.  Another 
body  of  work  has  been  concerned  with  the  effects  of  uplink 
noise  passing  through  a  nonlinearity  [14] -[24],  however,  the 
signal  is  presumed  to  be  wideband  (no  ISI)  and  the  non- 
linearities  are  treated  as  hard-  or  soft-limiters.  Some 
[25] -[26]  have  merely  developed  characterizations  of  typical 


nonlinearities  and  of  those  few  that  have  dealt  with  both 
ISI  and  uplink  noise  preceding  a  nonlinear  transponder 
[27]- [30],  most  have  used  just  some  form  of  limiter  and  not 
a  more  realistic  transponder  model. 

With  the  growth  of  the  satellite  field  bringing  more 
stringent  design  specifications,  there  has  been  a  realization 
that  the  flexibility  inherent  in  Monte  Carlo  simulation  is 
more  important  than  the  inherent  speed  of  the  analytical 
methods.  Coupled  with  the  ability  of  simulation  to  achieve 
any  level  of  accuracy  desired,  this  technique  is  seen  to  be 
a  more  powerful  and  desirable  one  than  performance  prediction 
using  purely  analytical  methods.  Since  it  has  been  shown 
[31]- [37]  that  Monte  Carlo  techniques  can  be  used  to  deal  with 
intersymbol  interference  in  the  absence  of  uplink  noise,  the 
prediction  via  simulation  problem  boils  down  to  that  of  how 
to  consider  the  effects  of  the  uplink  noise  without  taking  an 
unacceptable  amount  of  computational  time  to  do  so.  Modification 
to  the  noise  statistics  via  importance  sampling  [43] -[46]  is 
one  approach  that  has  been  taken  to  reduce  the  amount  of  run 
time.  Another  involves  the  estimation  of  the  tail  probabilities 
of  the  output  signal  level  using  either  extreme  value  distri¬ 
bution  theories  or  those  of  asymptotic  approximation  [38] -[42], 

As  far  as  the  hybrid  technique  goes  [2],  it  merely 
consists  of  using  a  mathematical  model  for  the  nonlinearity 
while  retaining  the  Monte  Carlo  model  for  the  rest  of  the 
system.  This  has  the  flexibility  problem  of  any  analytical 
method  in  that  a  highly  accurate  model  of  the  nonlinearity 
eliminates  the  ability  to  compare  differing  transponder 
configurations  and  components. 

The  remainder  of  this  thesis  is  divided  into  four  parts. 

In  the  next  section,  a  review  of  the  ideas  that  form  the 
basis  of  system  error  performance  prediction  and  its  importance 
is  presented.  Chapter  3  examines  the  status  of  performance 
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prediction  through  a  discussion  of  the  methods  mentioned 
in  the  preceding  paragraphs.  An  analytical  examination  of 
a  potential  prediction  tool’s  effectiveness  is  conducted  in 
Chapter  4  and  the  implications  and  conclusions  of  the  study 
are  presented  in  the  final  section. 


• » > 


STATEMENT  OF  THE  PROBLEM 


2.0 


A  typical  satellite  communications  channel  is  a  bandpass 
system  with  a  center  frequency  in  the  4-20  GHz  range  and  a 
bandwidth  of  500  MHz.  For  notational  convenience  in  dis¬ 
cussing  this  system,  we  shall  make  use  of  the  well-known 
fact  that  a  bandpass  system  can  be  modeled  by  an  equivalent 
baseband  system  through  utilization  of  the  concept  of  the 
complex  envelope  representations  of  bandpass  signals.  This 
modeling  is  effective,  of  course,  only  if  the  performance  of 
the  system  does  not  exhibit  frequency  dependence  above  and 
beyond  that  of  atmospheric  attenuation,  fading,  etc.  Since 
the  prediction  of  the  bit  error  probability  in  digital 
communications  is  almost  entirely  dependent  on  determining 
the  effects  of  noise  and  intersymbol  interference,  both  of 
which  are  easily  modeled  in  the  baseband  with  no  loss  of 
generality,  we  are  free  to  use  a  baseband  equivalent  system 
in  all  future  discussions.  Consequently,  Figure  1  presents 
a  generalized  model  of  the  typical  satellite  communications 
system  of  this  study.  It  should  be  noticed  that  the  notation 
used  in  Figure  1  and  in  the  following  paragraphs  adheres 
closely  to  that  of  Shanmugam  and  Balaban  in  [45]. 


2.1  The  System 

The  system  being  modeled  is  seen  to  represent  a  dedicated, 
single  carrier  per  transponder  configuration  that  can  easily 
be  thought  of  as  simulating  time-division  multiple-access 
(TDMA)  techniques.  Of  the  five  functional  areas  into  which 
the  system  is  grouped,  the  first  to  be  seen  by  a  user  is  the 


Typical  Satellite  Conmtmications  System 
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transmitter,  consisting  of  the  modulator  and  a  bandlimiting 
filter.  The  filter  is  a  conglomerate  representation  of  the 
various  stages  of  filtering  encountered  in  an  actual  transmitter, 
such  as  those  for  pulse-shaping  (a  function  that  is  performed 
entirely  in  the  transmitter  in  most  satellite  systems)  and 
for  fitting  the  signal  to  the  FCC  transmission  mask. 
"Bandlimiting"  simply  indicates  that  the  filter  has  the 
narrowest  bandwidth  (and  therefore,  smallest  bandwidth-signal 
pulse  time  (BT)  product)  of  any  in  the  system  and  that  all 
significant  xSI  arises  at  this  filter. 

Following  the  transmitter  is  the  section  of  the  system 
known  as  the  uplink.  This  is  a  representation  of  the  effects 
of  transmitting  a  signal  through  the  earth's  atmosphere  to 
the  communications  satellite.  Under  normal  operating  conditions, 
this  functional  area  within  the  system  may  be  characterized  as 
a  linear  channel,  or  filter,  with  additive  Gaussian  noise  at 
the  output.  Since  the  downlink  part  of  the  system  is  a  duplicate 
of  the  uplink,  it  can  be  characterized  in  the  same  manner. 
Sandwiched  between  the  earth-satellite  links  is  that  group 
of  components  that  are  comprised  by  the  satellite  transponder. 
This  group  consists  of  two,  almost  identical  filters,  one 
acting  as  the  transponder  receive  filter  with  the  other 
filling  the  role  of  the  transmit  filter,  and  a  travelling- 
wave  tube  amplifier.  Due  to  the  wide  variety  of  modems  that 
are  used  in  the  earth  stations  (i.e.,  the  lack  of  standardi¬ 
zation  of  telecommunications/data  communications  equipment) , 
the  filters  are  required  to  perform  their  duties  yet  be 
compatible  with  most  earth  stations.  In  addition,  the  unique 
circumstances  entailed  by  satellite  operation  require  that  the 
filters  be  very  dependable  while  economic  factors  require 
minimization  of  the  cost  of  the  system.  Hence,  the  filters 
are  as  nearly  identical  as  possible,  relatively  wideband  in 
nature,  and  perform  their  functions — restricting  the  effective 
bandwidth  of  the  uplink  noise  for  one  and  eliminating  any 
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out-of-band  signals  created  by  the  TWTA  for  the  other— 
without  introducing  any  time  dispersion  or  other  significant 
signal  degradation.  It  is  at  this  point  in  the  system  that 
all  of  the  nonlinear  operation  is  introduced.  Although  the 
high  power  amplifier  (HPA)  used  to  boost  the  level  of  the 
transmitted  signal  is  frequently  operated  in  its  nonlinear 
region,  the  effect  of  the  TWTA  is  dominant  and,  in  any  case, 
suffices  to  address  the  problem  of  predicting  the  performance 
of  a  typical  satellite  system. 

In  Figure  2  the  input-output  curves  for  a  typical  TWTA 
are  shown.  Although  there  is  a  linear  region  of  operation, 
it  is  easy  to  see  that  maximum  power  output,  and  thereby 
maximum  downlink  Eb/N0  (energy  per  bit/noise  power  spectral 
density) ,  is  achieved  when  operation  is  in  the  nonlinear 
region.  For  this  study,  the  TWTA  shall  be  considered  to 
operate  always  at  saturation  as  is  usually  the  case  in 
realistic  TDMA  systems.  Looking  past  the  transponder  one 
sees  the  downlink,  already  noted  as  being  identical  to  the 
uplink,  and  the  final  section  of  the  system,  the  receiver. 
This  section  is  a  mirror  image  of  the  transmitter  with  two 
important  exceptions.  First,  the  composite  receive  filter  is 
similar  in  function  to  the  transmit  filter  but  does  not 
bandlimit  the  signal  quite  so  much  therefore  introducing  less 
ISI,  and  second,  there  is  obviously  some  form  of  sampling 
and  decision  logic  integrated  with  the  demodulator.  Thus, 
the  receiver  input  is  filtered  to  restrict  its  bandwidth  and 
then  passed  through  the  demodulator/decision  logic  from  which 
the  information  reaches  its  destination. 


I 

I 

! 

i 


Now  suppose  that  one  wishes  to  estimate  the  probability 
of  error  for  the  system  just  described.  We  will  only  look 
at  the  case  of  binary  pulse  amplitude  modulation  (PAM)  since 
this  reduces  the  complex  envelope  to  a  real  signal  and  is 
sufficient  for  illustrative  purposes.  It  can  be  assumed  that 
the  input  to  the  system  is  a  sequence  of  binary  dat  {b^} 
having  a  bit  duration  of  T^.  Following  modulation,  the  input 
to  the  transmit  filter  has  the  general  form 


D ( t)  =  k  afreet (t-kTb) 


rect(t-kTb)  = 


1  kTfa  <  t  <  (k+l)Tb 
0  elsewhere 


where  the  rect(t)  pulse  was  arbitrarily  chosen  for  computational 
ease  and  {ak>  is  the  modulating  amplitude  sequence  with 
a^  =  -A  or  +A  when  b^  *  0  or  1  respectively.  With  the  transmit 
filter  assumed  to  have  the  transfer  function  ht(t)  and  the 
impulse  response 


Ht(t)  =  ht(t)cos  wQt, 


the  transmitter  output,  S(t),  can  be  expressed  as 


S(t)  =  k  akp(t-kTb) 


S(kTb)  =  akp{0)  +  ^  ak_mp[(k-m)Tb] 


where  p(t)  is  the  shaped,  bandlimited  signal  pulse 


p(t)  =  0.5[rect(t)  *  ht(t)]  (6) 

In  equation  (5) ,  the  first  term  on  the  right  represents  the 
k-th  transmitted  bit  while  the  summation  represents  the 
residual  effect  of  all  other  transmitted  bits.  This  second 
term  represents  the  ISI  that  is  caused  by  the  bandlimiting 
filter.  As  it  passes  through  the  uplink,  the  signal  S(t) 
is  corrupted  by  noise,  n(t),  that  can  be  modeled  by  a  zero- 
mean  Gaussian  random  process.  Representing  the  thermal  noise 
that  degrades  the  signal  in  the  uplink,  n(t)  is  additive  in 
nature,  and  hence  the  transponder  input,  X(t),  can  be  expressed 
by 


X  ( t)  =  S  (t)  +  n  ( t) 

=  £  akp(t-kTb)  +  n (t) 


(7) 


Passing  through  the  filtering  and  nonlinear  transformation 
of  the  transponder,  the  signal  as  it  enters  the  downlink  is 
described  by 


z(t)  =  f  [X  ( t)  j  (8) 

where  f(x)  is  the  input-output  relationship  of  the  transponder. 
It  is  beyond  this  point  that  the  probability  of  error  can  be 
seen  to  have  been  rendered  very  difficult  to  track  analytically 
since  the  statistics  of  both  the  uplink  noise  and  the  ISI  have 
been  altered.  In  fact,  the  alteration  is  frequently  such 
that  it  is  impossible  to  separate  the  effects  the  two  have 
on  the  information-bearing  part  of  the  signal. 


Within  the  downlink,  the  signal  is  again  corrupted  by 
a  random  process  representing  AWGN,  v(t),  so  that  the  input 
to  the  receive  filter,  U(t),'  has  the  form 

U(t)  *  Z(t)  f  v ( t )  .  (9) 

If  the  transfer  function  of  the  receive  filter  is  assumed  to 
be  hr(t)  with  the  impulse  response 

Hr(t)  =  hr(t)cos  ui0t,  (10) 


then 


Y(t)  -  0 . 5  [U  ( t)  *  hr(t)] 

=  0.5  ( [Z  (t)  *  hr(t)]  +  [v  ( t)  *  hr(t)])  (11) 

=  J  Z*  (t-nTb)  +  v*  (t) 

is  the  filtered  input  to  the  demodulator  and  decision  logic. 

It  is  easy  to  see  the  complications  in  predicting  the  system 
BER  if  one  notes  that  (11)  contains  ISI  that  comprises  the 
effects  of  a  nonlinear  transformation  of  the  desired  signal 
+  noise  +  transmit  filter  ISI.  For  reasons  that  shall  become 
obvious,  we  will  assume  that  the  statistical  properties  of 
the  system  inputs,  n(t)  and  v(t) ,  as  well  as  the  system  transfer 
characteristic — four  filter  transfer  functions  and  the  TWTA 
transfer  characteristic — are  known. 

At  this  stage,  one  is  confronted  with  the  actual  attempt 
to  estimate  the  BER,  or  P2(e),  for  the  system.  This  quantity 
is  given  by 


00 

p2(e)  =  (Jf  qY|bk=0(y)dy}PrCbks°] 
< 

*  <aYib,  =i<y)dy  pr[bk=i] 


(12) 


where  <  is  the  decision  threshold  (typically  0  for  the  system 
modeled  by  Figure  1)  and  qy^^Cy)  is  the  conditional  probability 
density  function  of  the  output  Y,  given  the  input  b^.  The 
performance  prediction  problem  is  thus  seen  to  be  that  of 
estimating  the  two  conditional  pdf's.  This  task  is  made 
quite  difficult  by  the  nonlinear  transformation  that  has 
altered  the  statistics  of  the  transmitted  signal  and  the  uplink 
noise.  Additionally,  an  analytical  approach  is  inherently 
inflexible;  there  is  no  room  for  "black-box"  modeling  of  the 
system  and  its  components;  and  is  rendered  somewhat  inaccurate 
by  the  need  to  approximate  the  effects  of  the  system  non¬ 
linearity  and  the  ISI.  On  the  other  hand,  Monte  Carlo 
simulation  estimates  the  pdf's  by  what  is  essentially  the 
creation  of  two  histograms.  This  allows  for  great  flexibility 
and  the  achievement  of  any  level  of  accuracy  desired — although 
the  effects  of  ISI  must  still  be  approximated  by  a  finite 
number  of  pulses.  There  is  a  price  that  is  paid  for  the 
accuracy,  however.  When  there  is  a  random  process  preceding 
the  TWTA,  such  as  uplink  noise,  enough  samples  must  be  examined 
to  allow  the  process's  statistics  full  reign;  this  number 
is  on  the  order  of  10/P(e)  [44].  Since  a  typical  BER  for  a 
satellite  system  is  one  error  out  of  every  10®  symbols,  Monte 
Carlo  methods  require  about  10?  symbols  and  hence,  an 
unacceptably  large  amount  of  computational  time. 

Obviously,  the  main  difficulty  in  the  estimation  of  the 
probability  of  error  arises  due  to  the  fact  that  intersymbol 
interference  and,  more  importantly,  noise  occur  prior  to  the 
nonlinear  amplifier.  Just  what  effects  this  has  and  why  it 
is  difficult,  and  necessary,  to  handle  will  be  discussed  in 
the  next  section. 
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2. 3  A  Discussion  of  the  Problem 

In  an  ideal  digital  communications  system,  the  pulse¬ 
shaping  filters  are  assumed  to  be  ideally  matched  thus 
eliminating  any  possibility  of  the  occurrence  of  intersymbol 
interference.  Even  when  the  filters  cannot  be  matched,  as 
in  the  typical  satellite  communications  case  where  lack  of 
standardization  forces  all  pulse-shaping  to  be  done  at  the 
transmitter,  the  vast  majority  of  the  ISI  can  be  eliminated 
through  the  usage  of  a  good  equalizer— if  the  channel  and 
filter  characteristics  are  known  exactly.  Since  that  condition 
is  impossible  to  meet,  there  will  always  be  some  ISI  in  a 
real,  bandlimited  system;  there  has,  in  fact,  been  a  great 
amount  of  work  done  exclusively  on  the  problem  of  determining 
the  effect  of  the  ISI  on  a  system's  performance  [3]— [15] . 

Most  of  the  work  already  done,  however,  has  dealt  only  with 
linear  channels  (with  additive  Gaussian  noise,  of  course) 
and  without  this  restriction  the  solution  to  the  problem  has 
not  been  easily  found. 

ISI,  the  time-dispersion  of  a  signal  pulse's  energy  that 
causes  the  amplitude  of  a  signal  at  any  given  moment  to  be 
dependent  on  several  different  pulses,  combines  with  thermal 
noise  to  be  the  primary  causes  of  errors  in  a  digital  system. 

In  a  linear  channel  with  linear  filters,  the  ISI  has  an  effect 
similar  to  that  of  noise  but  non-Gaussian  in  nature.  The 
exact  analysis  of  the  joint  effects  of  noise  and  intersymbol 
interference  involves  the  computation  of  a  very  complex 
probability  distribution  [4], [5].  For  practical  systems,  in 
fact,  an  exact  analysis  can  seldom  be  carried  out  [6]. 

With  a  little  thought,  it  can  be  seen  that  the  nonlinear 
characteristic  of  the  typical  satellite  channel  changes  the 
"seldom"  of  the  preceding  sentence  to  "probably  never."  A 
nonlinear  transformation  of  the  uplink  ISI  changes  its  statistics. 
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As  a  result,  only  some  of  the  IS I  can  be  removed  by  the  matched 
filter  at  the  receiver.  If  we  assume  that  ISI  arises  on  the 
downlink  as  well,  then  the  analytical  problem  becomes  even 
more  of  a  mess  and  attempts  at  finding  anything  other  than 
bounds  quickly  become  hopeless. 

There  are  ways  to  approximate  the  effects  of  intersymbol 
interference.  Analytically,  one  can  find  worst-case,  best-case 
bounds,  or,  using  various  series  expansions,  one  can  find 
somewhat  tighter  probability  of  error  bounds.  These  are, 
however,  still  bounds  and  not  the  fairly  accurate  estimations 
that  are  desired.  A  way  does  exist  to  achieve  a  fairly 
accurate  estimate  of  the  effects  of  ISI;  The  enumeration  of 
all  possible  values  the  interference  can  have  and  the  subsequent 
averaging  of  the  P(e)  for  each  value.  While  this  can  be  done 
analytically,  the  enumeration  method  tends  to  be  very  cumbersome 
in  most  cases  (even  with  the  finite  number  of  interfering 
pulses  that  must  be  assumed) ,  and  is  much  better  suited  to 
simulation.  Studies  have  been  made  that  show  enumeration  in 
conjunction  with  Monte  Carlo  simulation  to  be  quite  adequate 
in  estimating  the  effects  of  ISI  on  system  performance. 

Since  the  argument  has  already  been  made  convincing  us 
of  the  desirability  of  simulation  over  the  analytical  methods 
as  a  performance  prediction  tool ,  we  can  consider  the  question 
of  how  one  determines  the  effects  of  ISI  on  the  system  BER 
to  be  answered  via  enumeration.  That  leaves  only  the  uplink 
noise  as  a  major  stumbling  block  in  specifying  a  prediction 
tool  for  use  in  designing  satellite  systems.  In  many  of  the 
early  communications  satellites,  enough  power  was  broadcast 
from  the  transmitter  to  make  the  uplink  E^/Nq  significantly 
larger  than  that  of  the  downlink.  This  situation  meant  that 
only  a  very  insignificant  number  of  errors  were  caused  by  the 
uplink  noise.  Thus,  the  methods  that  have  been  developed  to 
estimate  the  BER  have  been  able  to  assume  that  all  of  the 


17 


thermal  noise  affecting  a  system  arises  at  the  input  to  the 
receiver.  Alternately,  if  the  TWTA  could  be  considered  to 
be  operating  in  its  linear  region  to  avoid  intermodulation 
(as  was  the  case  for  most  of  the  early,  multi-carrier, 
frequency-division  multiple  access  (FDMA)  satellites) ,  the 
uplink  noise  was  combined  with  the  downlink  noise  either 
linearly  or  in  some  other  mathematically  simple  manner.  This 
approach  tended  to  yield  the  same  effective  result  as  did 
ignoring  the  uplink  noise  entirely. 

In  the  past  few  years,  however,  the  number  of  satellites 
and  users  has  grown  dramatically  while  the  usable  spectrum 
has  not.  More  and  more,  the  multiple  access  technique  of  choice 
and  necessity  is  TDMA.  This  allows  the  transponder  to  be 
operated  at  a  much  more  power-efficient,  albeit  nonlinear, 
level.  Operating  the  TWTA  in  its  nonlinear  range  means  that, 
if  the  uplink  noise  is  significant  at  all,  the  (Eb/NQ)U 
cannot  be  directly  combined  with  the  (Eb/N0)d  with  any  hope 
for  accuracy — and  the  uplink  noise  is  frequently  significant. 

The  necessity  of  using  TDMA  now  is  that  the  increased  number 
of  satellites,  and  hence,  proximity,  require  that  less  power 
be  broadcast  to  each  satellite  to  avoid  interfering  with 
transmissions  to  nearby  satellites.  Additionally,  there  has 
been  a  virtual  explosion  in  the  number  of  small  earth  stations 
that  are  being  used,  thus  compounding  the  interference  problem 
and  requiring  an  even  larger  reduction  in  the  transmitted 
power.  This  reduction  in  transmit  power  has  been  great 
enough  that  (Eb/NQ)U  is  now  of  roughly  the  same  order  as 
that  for  the  downlink  thereby  necessitating  its  consideration 
as  a  major  cause  of  error. 

The  nonlinear  transformation  does  more  than  just  preclude 
the  direct  combination  of  uplink  and  downlink  Eb/NQ's.  It 
forces  one  who  is  taking  an  analytical  approach  to  the  problem 
to  calculate  the  nonlinearly  transformed  pdf  of  the  noise  at 


the  output  of  the  TWTA;  aside  from  being  a  rather  difficult 
undertaking  in  and  of  itself,  it  is  easy  to  see  that  the  pdf 
will  be  heavily  dependent  upon  the  characteristic  of  the 
particular  nonlinearity  thereby  placing  severe  limitations  on 
the  flexibility  of  the  method  used  as  a  performance  prediction 
tool.  In  addition,  most  of  the  work  that  has  been  done  on 
this  problem  has  either  ignored  the  existence  of  intersymbol 
interference  [14] -[26]  and/or  has  treated  the  TWTA  nonlinearity 
as  some  kind  of  limiter  [27] -[30]. 

It  is  obvious  that  the  IS I  and  the  uplink  noise  must  both 
be  considered  if  one  wishes  to  derive  an  accurate  estimate  of 
a  system's  BER.  Passage  through  the  TWTA  causes  the  two  to 
"freeze"  one  another  into  the  signal  in  such  a  way  that  the 
individual  effects  cannot  be  separated  from  each  other  or  from 
the  information-bearing  part  of  the  signal.  The  difficulties 
of  using  analytical  methods  on  either  of  the  effects  standing 
alone  has  already  been  stated.  Together,  the  inevitable 
conclusion  is  that  the  analytical  methods  will  be  inadequate. 

Monte  Carlo  simulation  encounters  its  own  problems  when 
dealing  with  the  uplink  noise,  although  handling  the  ISI  as 
well  as  the  uplink  noise  is  not  one  of  them.  As  mentioned 
previously,  a  typical  satellite  system  BER  is  on  the  order  of 
10”6  meaning  that  least  10 7  bits  need  to  be  run  through  the 
simulation  to  give  an  accurate  error  count.  The  computational 
time  involved  in  handling  this  number  of  bits  makes  regular 
Monte  Carlo  simulation  unacceptable  as  a  performance  prediction 
tool. 

It  immediately  becomes  obvious  that  there  is  a  need  for 
some  modification  of  the  basic  Monte  Carlo  methods  that  will 
greatly  reduce  the  number  of  samples  required  to  acquire  a 
good  estimate  of  the  BER  of  a  communications  system.  In  the 
following  chapters,  a  look  at  many  of  the  previous  attempts 
at  developing  some  form  of  performance  prediction  tool 
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(analytical  or  simulation)  will  be  presented  followed  by  a 
study  of  two  techniques  that  appear  to  hold  the  greatest 
promise  of  solving  the  problem.  In  addition,  a  combination 
of  these  two  techniques  will  be  studied  to  determine  what  is 
the  best  performance  prediction  tool  existing  at  this  time. 


SURVEY  OF  THE  PREVIOUS  WORK 


3.0 


From  the  previous  discussions  we  know  that  the  main 
causes  of  detection  error  in  a  digital  communications  system 
are  intersymbol  interference  and  additive  noise.  When  a 
system  is  afflicted  by  the  two  of  these  together,  exact 
estimation  of  the  error  probability  becomes  computationally 
complex.  If,  further,  a  typical  satellite  system  such  as 
that  shown  in  Figure  1  is  considered,  one  is  faced  with  the 
situation  where  the  transmitted  signal  is  corrupted  by  ISI 
and  noise  prior  to  passing  through  a  TWTA  operating  in  its 
nonlinear  region.  In  addition,  before  the  decision  logic  is 
reached,  the  signal  is  corrupted  by  even  more  ISI  and  additional 
noise.  Under  these  circumstances,  it  is  easy  to  imagine  that 
the  already  quite  complex  task  of  estimating  the  system  BER 
becomes  almost  impossible. 

The  necessity  of  having  a  performance  prediction  tool, 
however,  is  such  that  a  large  number  of  approaches  towards 
solving  the  error  estimation  problem  have  been  developed.  It 
has  been  mentioned  previously  that  there  are,  primarily,  two 
categories  into  which  the  various  prediction  tools  may  be 
classified.  These  categories  are  the  result  of  the  fact  that 
there  are  two  unreconcilable  requirements  that  an  ideal 
performance  predictor  should  meet.  The  first  of  the  require¬ 
ments  is  that  the  prediction  tool  be  very  flexible  with 
regard  to  the  system  components,  configuration,  modulation 
formats,  etc.  This  has  led  to  the  growth  of  Monte  Carlo 
simulation,  in  which  the  system  is  just  a  set  of  "black  boxes" 
that  can  be  arranged  in  any  desired  manner  and,  further,  the 
signals  that  enter  the  system  are  treated  as  random  processes. 
The  second  requirement  is  that  the  performance  prediction  tool 


deliver  answers  speedily  and  with  a  great  deal  of  accuracy. 

An  obvious  outgrowth  of  this  requirement  has  been  the 
development  of  various  analytical  methods  that  attempt  to 
develop  a  unified  deterministic  system  model  with  a  set  of 
parameters  from  which  the  performance  can  be  calculated. 

While  there  are  some  hybrid  methods  that  employ  character¬ 
istics  of  each  of  the  aforementioned  techniques,  the  vast 
majority,  if  not  all,  may  be  placed  into  either  the  simulation 
category  or  the  analytical  methods  category.  With  this  in 
mind,  the  following  sections  will  provide  a  description  of 
the  fundamentals  of  each  different  type  of  approach  taken  in 
the  attempt  to  find  an  "ideal"  performance  prediction  tool. 

A  synopsis  of  the  advantages  and  disadvantages  accruing  from 
each  approach  shall  also  be  presented.  Most  of  the  literature 
reviewed,  it  will  probably  be  noted,  is  specifically  applicable 
to  some  type  of  phase-shift-keyed,  (PSK) ,  system.  Done 
primarily  because  of  the  pervasiveness  of  PSK  as  a  modulation 
format,  there  is  no  loss  of  generality  in  the  analyses 
presented  as,  on  careful  consideration,  one  realizes  that 
almost  any  modulation  format  can  be  expressed  via  PSK.  In 
the  following,  the  use  of  PSK  as  the  operational  modulation 
of  discussion  will  be  adopted. 


3.1 


The  Analytical  Methods 


In  Figure  3,  a  modified  version  of  the  digital  satellite 
communications  system  model  shown  in  Figure  1  is  presented 
for  use  in  discussing  the  various  analytical  approaches 
taken  towards  solving  the  problem  of  accurately  predicting 
a  system's  BER.  Although  the  two  models  are  essentially  the 
same,  the  signal  representations  of  Figure  1  have  been  changed 
to  accommodate  the  fact  that  Figure  3  is  meant  to  represent 
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a  bandpass  PSK  system  similar  to  those  discussed  in  the 
literature.  Consequently,  a  brief  run-through  of  the  system's 
operation  is  in  order. 

An  M-ary  PSK  system  operates  by  modulating  its  carrier 
with  the  phase  sequence  {$j_}  where  each  transmitted  phase 
is  assumed  to  take  one  of  M  equally  probable  values,  typically 

*k  =  M  (1+2 A)  l  =  0,1,..., M-l  (M  even)  (13) 

If  one  assumes  that  the  input  to  the  system  is  binary  data, 

then  each  phase  angle  <j>k  cam  be  considered  as  representing  an 

m  bit  data  symbol  where  M  =  2m.  Accordingly,  the  symbol 

duration  is  T_  with  T  =  mT.  where  T.  is  the  bit  duration, 
s  s  b  b 

Keeping  this  information  in  mind,  we  shall  proceed  to 
describe  and  follow  an  M-ary  PSK  signal  as  it  passes  through 
the  system  of  Figure  3.  A  phase  modulated  signal  may  be 
thought  of  in  one  of  two  ways:  Either  as  a  constant  amplitude 
carrier  to  which  the  sequence  { >  has  been  added,  or,  as 
two  amplitude  modulated  carriers  with  one  of  the  carriers 
delayed  90°  relative  to  the  other.  In  the  former  case,  the 
signal  at  the  output  of  the  modvilator,  Dp(t),  would  be  written 

Dp  ( t)  =  A  cos  [u>Qt  +  a  (t)  ]  (14) 

where  A  is  an  arbitrary  constant,  is  the  center  frequency 
of  the  carrier  (in  radians/Hz)  and  a(t)  is  the  modulating 
phase  angle.  a(t)  can  be  represented  as 

ct  ( t )  =*  £  $m(t-kTs)  rect  (t-kTg)  (15) 

with  the  rect(t-kT  )  pulse  defined  by  (2)  and  serving  merely 
for  reasons  of  computational  convenience  and  not  as  a  requisite. 
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As  a  matter  of  fact,  any  pulse  of  width  T  is  permissible. 
Bandlimiting,  however,  creates  a  cumbersome,  if  not  overly 
complex,  signal  representation  that  is  not  well-suited  to 
the  purposes  of  this  discussion.  We  will  therefore  favor  the 
latter  method  for  representing  PSK  signals.  In  this  case, 

Dp(t)  =  £  rect ( t-kTs)  cos  <fr{t-kTg)  cos  uQt  -  rect(t-kTs) 

sin  $  (t-kT  )  sin  o»_t  (16) 

3  O 

=  £  rect(t-kTg)  cos  [a>0t  +  <j>  (t-kTg)] 

and  once  again,  rect(t-kTg)  is  chosen  for  convenience's  sake. 
Throughout  the  rest  of  this  chapter,  please  note  that  the 
time  dependence  of  all  variables  will  be  ignored  whenever 
possible  once  they  have  been  introduced  and  defined. 

The  modulated  signal  now  passes  through  the  bandlimiting 
transmit  filter  that  is  assumed  to  have  the  impulse  response 
given  by  (3).  From  this  filter,  the  transmitter  output,  S^Ct), 
emerges  with  the  form 


sp(t) 

*  £  p(t-kTg)  cos  [u0t  +  <j>  (t-kTg)  ] 

(17)  ; 

p(t) 

*  H  [rect (t)  *  hfc (t) ] 

(18) 

and  for  kT  <.  t 
s 

<  (k+1)  T 

5 

Sp  *  P(t-kTs) 

cos  [uiQt  +  $  ( t-kTg )]  +  jJkP<t-iTg) 

cos  [ aiQ t+<|>  (t-iTg)] 

=  (Pj^COS$j^ 

+  i^kPi  cos*.)  cos  u)Qt  -  <Pk  sin 

♦k  +  i5k*l  sin  *i> 

sin  !iiQt 

(19) 

=  A (t-kTg) 

cos  [u_,t  +  Ht-kT  )] 

o  s 
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which  means  that 


A(t-kTs) 


=  A,. 


=  [  (pk  cos  $k  +  cos 

+  (pk  sin  $k  +  sin 


(20) 


Ip  (  t-KT^ 


* 


=  arctan 


pk  sin<jik  +  i7ikPi  sin*. 
pk  cos*k  +  iJkPi  sin*J 


(21) 


where  Pk  =  p(t-kTg)  and  p^  =  p(t-iTg).  If  the  time  assumption, 
kT  <.  t  <  (k+l)T  is  dropped,  then  S  may  again  be  expressed 
as 


Sp  =  k  Ak  cos^ot  +  V 


(22) 


where  it  is  seen  that  Ak  =  pk  and  *k  =  <t>k .  In  equations 
(19)-(22),  the  ISI  effects  are  shown  explicitly  and  then 
absorbed  back  into  the  simplified  expression  of  (17)  . 

while  traversing  the  uplink  section  of  the  channel,  the 
signal,  Sp,  is  corrupted  by  noise,  n(t),  having  the  form 


n(t)  =  n  (t)  cos  w  t  -  n  (t)  sin  w  t. 

G  OS  O 


(23) 


In  the  preceding  equation,  nc  and  ns ,  the  inphase  and  quadra¬ 
ture  components  of  the  noise,  are  independent,  stationary, 


zero-mean  Gaussian  random  processes  with  each  having  a  variance 
on^.  Since  this  uplink  noise  happens  to  be  additive  in 
nature,  the  signal  that  enters  the  transponder,  Xp(t),  may 
be  expressed  by 


xp(t> 


S  +  n 
P 


k  k 


A.  cos(a)  t  +  (JO  +  n_  cos  to  t  - 


n_  sm  to  t 
s  o 


k  Rk(t)  cos  [ (oQ t  +  Sk  (t)  ] 


Rk(t)  =  [  (Ak  cos  ijjk  +  nc) 


+  {Ak  sin  ^k  +  rig)2]15 


6k(t) 


sin  i(/k  +  n£ 


Ak  Sln  ^k 
arCta"  cos 


with  Bk  obviously  representing  the  phase  effects  of  the 
uplink  noise  and  ISI,  as  well  as  the  transmitted  phase. 

From  the  set  of  equations  presented  above,  along  with 
( 19 )  —  ( 22 ) ,  it  is  easy  to  see  that  the  original  PSK  signal 
has  been  greatly  altered.  The  overall  phase  and  amplitude 
have  become  heavily  dependent  upon  many  different  random 
processes,  meaning  that  the  signal  presented  to  the  transponder 
has  a  variable  amplitude  across  each  symbol — a  situation  that 
can  have  much  significance  vis*  PM(e),  the  error  probability. 

As  mentioned  earlier,  the  typical  TDMA  satellite  communi¬ 
cations  system  operates  with  its  transponder  delivering 
maximum  possible  output  power.  This,  of  course,  means  that 
the  TWT  amplifier  forming  the  heart  of  the  transponder  is 


operated  within  its  nonlinear  region  (recall  Figure  2) .  The 
resulting  AM-AM  and  AM-PM  conversion  effects,  expressed  as 
f(*)  and  g(*)  respectively,  mean  that  the  varying  amplitude 
of  the  input  signal  produces  a  randomly  variable  phase  error 
rather  than  the  predictable  one  of  an  ideal  PSK  signal. 
Consequently,  the  IS I  and  uplink  noise  find  that  they  have 
"frozen”  one  another  into  the  signal.  The  signal  that  is 
output  to  the  downlink  from  the  transponder,  Z  (t) ,  can 
therefore  be  described  by 

Zp(t)  =  k  f(Rk}  c°s[u>0t  +  8k  -  g(iy]  (27) 

It  is,  at  this  point,  already  easy  to  comprehend  that  the 
effects  that  the  uplink  noise  and  IS I  have  on  the  system  BER 
has  severely  inhibited  its  analytical  tractability .  Not  only 
have  the  statistics  of  those  two  primary  causes  of  signal 
degradation  been  altered  nonlinear ly,  but  they  have  become 
increasingly  intertwined  with  one  another  and  with  the  desired 
part  of  the  signal.  So  much  so,  in  fact,  that  it  is  impossible 
to  assess  the  impact  of  one  or  the  other  by  itself. 

As  if  the  analytical  problem  hasn't  been  made  difficult 
enough  at  this  point,  the  signal  is  once  again  corrupted  by 
noise.  Like  that  of  the  uplink,  the  downlink  noise,  v(t) , 
has  the  form 

v(t)  =  v_(t)  cos  U)  t  -  v_(t)  sin  u^t  (28) 

c  os  o 

where  v  and  v  are  equivalent  to  their  uplink  counterparts, 

c  s  2 
nc  and  ng,  although  having  a  unique  variance  av  which  may  or 

may  not  be  the  same  as  a  Since  v  is  also  additive  in 

nature,  the  input  to  the  receiver,  U  (t)  ,  can  be  written  as 


Up(t)  =  zp  +  V 


*  k  cos[uot  +  8k  -  g(Rk)]  +  vc  cos  uiQt 


-  vc  sm  i^t 
s  o 


k  "k 


Wk(t)  COs[(JQt  +  0k(t)] 


(t)  =  [(f(\)  cos[8k-g(Rk)  ]  +  vc) 2 

+  (f  (ly  sint  8k-g(Rk)  ]  +  vs)2Jis 


‘f  0^)  sin  8k 


9k(t)  =  arctan  ^ 


Hi] 

+  VnJ 


with  9k  representing  the  cumulative  phase  at  the  receiver 
input. 

Passing  into  the  receiver,  the  signal  once  again  encounters 
a  bandlimiting  filter.  Although  not  quite  as  narrow  as  the 
bandwidth  of  the  transmit  filter,  the  receive  filter  bandwidth 
is  still  sufficiently  small  to  introduce  some  ISI.  With  the 
assumption  that  this  filter  has  the  impulse  response  given  by 
(10),  the  filter  output,  Y  ( t ) ,  has  the  form 


YP(t)  =  k  j  Pk^Tg)  C°S(u)0t  +  9k) 


Pk(t)  =  Jj[wk(t)  *  hr(t)] 


Recalling  (19) ,  if  we  now  restrict  our  discussion  to  the 
time  period  jT  <  t  <  (j+l)Te,  the  filter  output  may  also 

5  S 

be  written 


*  “  »  *  •  *  •  •  •  *  »  “  *  *  »  “  k  '  »  *  *«’•*■*  * 

O  *  '  •••■**** 


Yp  =  k  pk(t_^Ts)  cos(«ot+0k)  +  k  nf/j  Pk(t-nTs)  cos(a,0t+6k)  (34) 

where,  as  usual,  the  second  term  on  the  right  side  represents 

the  ISI.  It  can  be  noted  that  Yp  is  a  double  summation  due 

to  the  ISI  introduced  by  the  receive  filter.  This  shows  the 

additive  nature  of  multiple  sources  of  ISI  in  series.  Yp  passes 

now  into  the  demodulator  where  it  is  demodulated  and  sampled 

and  where  the  decision  is  made  as  to  the  value  of  <*>y 

It  is  at  this  point  that  one  is  faced  with  the  problem 

of  analytically  determining  the  probability  of  error  for  the 

system.  If  we  denote  the  probability  density  function  of  Yp, 

given  that  the  transmitted  phase  <j>  =  as  q  ( Y  |  <(>  *  $.),  and 
£  r  J  3 

if  Oj  =  k  Pkft-jTg),  pn  “  k  pk^t“nTs^'  then  the  expression  for 

the  error  probability  is 

PM(e)  =  1  -  PrUj  -  gj(l+2j)  <  +  A.(t)  <  +  ^j(l+2j>] 


=  Pr [ | Aj (t) |  >  2m ( 1+2 j ) ] . 


Aj  (t) 


■  p- 

L 


sln  9k  +  n/i  °n 
cos  9k  +  n*j  pn 


sin  ek 
sin  0k 


Even  when  one  assumes,  as  we  do,  that  the  statistical  proper¬ 
ties  of  both  noise  functions  and  are  known,  the  determination 
of  the  pdf  of  A,  q(A),  is  an  extremely  complicated  process. 

The  development  of  the  various  analytical  methods  has 
followed  a  fairly  straightforward,  more  or  less  chronological 
path.  In  much  of  the  early  work,  even  before  satellite  systems 
became  the  primary  genre  under  consideration,  the  methods 
concentrated  on  how  to  determine  the  effects  of  ISI  in  con¬ 
junction  with  additive  noise.  The  systems  dealt  with  were 


primarily  linear  PAM  systems  with  additive  Gaussian  noise 
and  ISI  arising  either  because  of  mismatched  transmit  and 
receive  filters  or  other  reasons.  As  satellites  became  more 
common,  designers  were  faced  with  two-link  systems  containing 
a  nonlinearity.  Consequently,  the  next  step  along  the  path 
taken  considered  nonlinear  systems  with  additive  noise  on  both 
the  uplink  and  the  downlink.  These  approaches,  however, 
tended  to  ignore  the  existence  of  ISI.  Concurrently,  some 
researchers  examined  nonlinear  systems  with  intersymbol 
interference  but  no  noise  on  the  uplink.  Logically,  therefore, 
the  final  path  taken  has  been  the  consideration  of  systems 
containing  all  three  sources  of  signal  degradation.  That  is, 
the  consideration  of  systems  similar  in  nature  to  that  shown 
in  Figure  3  with  additive  noise  on  both  the  uplink  and  the 
downlink  and  a  nonlinear  satellite  transponder.  In  the 
following,  the  summaries  of  the  literature  will  be  grouped 
by  the  complexity  of  the  class  of  systems  with  which  they 
have  dealt. 

3.1.1  Approaches  to  Linear  Systems  with 

ISI  and  Additive  Noise 

At  the  lowest  level  of  complexity  are  those  approaches 
whose  sole  purpose  is  the  determination  of  the  effects  of 
intersymbol  interference  on  the  performance  of  linear  systems 
in  the  presence  of  additive  noise.  As  can  be  surmised,  the 
majority  of  this  work  was  done  at  a  time  when  FDMA  systems, 
with  their  resulting  linear  operation,  were  the  only  type 
being  used.  Nonetheless,  the  work  that  is  summarized  herein 
forms  the  basis  for  much  of  the  later  work  and  is,  therefore, 
important. 


Figure  4  shows  a  model  that  may  be  used  to  represent  a 
linear  system  with  additive  noise  and  ISI.  Given  that  this 
system  has  the  same  input  to  the  transmit  filter  as  the 

system  of  Figure  3,  then  the  receive  filter  output,  r(t), 

will  be  similar  in  form  to  that  of  Xp  described  previously. 
Namely,  for  kT  <  t  <  (k+l)T_,  and  with 

S  "  s 

Dp  =  ?  rect  (t-iTg)  cos(aiQt  +  ^)  ,  (14) 

r  (t)  =  p(t-kTg)  cos  U0t+<j>k)  +  ^^pft-iTg)  cos  (u^t+c^)  +  n 

=  [pk  cos  +  i^k  pi  cos  <t>±+  nc]  COS  a)Qt 

-  [pk  sin  <J>k  +  i?£k  p±  sin  4>i  +  nc]  sin  uQt 

=  r.  (t)  cos  a)  t  -  r_(t)  sin  t  (37) 

i  0  2  o 

where  n  is  the  same  as  that  of  (23)  and  p(t)  is 

p(t)  =  (rect (t)  *  [ht (t)  *  hc(t)  *  hr (t)  ] )  .  (38) 

In  (38)  above,  h  (t)  is  the  representation  of  the  effects 

C 

the  linear  channel  between  the  transmitter  and  the  receiver 
has  on  the  transmitted  signal. 

A  fundamental  approach  to  determining  the  effects  of  ISI 
in  any  type  of  system  is  that  of  finding  the  worst  case  bound 
An  upper  bound  to  the  probability  of  error  may  be  found  by 
assuming  that  the  ISI  always  takes  its  largest  value.  This 
bound  is,  in  general,  a  very  pessimistic  one  since  the  ISI 
only  assumes  its  largest  value  with  a  very  small  probability. 


Another  very  basic  method  of  determining  the  effects 
of  ISI  is  to  obtain  a  lower  bound  to  the  BER  by  considering 
only  a  finite  number  of  the  interfering  pulses.  Jones  [3] 
discusses  how  this  is  done  for  a  linear  PSK  system  containing 
ISI  and  AWGN.  His  conclucions  are  that  this  is  a  very 
cumbersome  method  of  calculating  ISI  effects  even  when  one  is 
only  dealing  with  the  effects  of  an  adjacent  pulse.  If, 
however,  the  ISI  due  to  distant  pulses  is  large,  one  comes 
upon  a  situation  wherein  the  bound  obtained  is  very  loose 
and  cannot  be  made  significantly  tighter  due  to  the  length 
of  the  required  computation. 

With  crude  upper  and  lower  bounds  realizable,  using 
worst  case  analysis  for  the  former  and  direct  enumeration  for 
the  latter,  the  obvious  next  step  down  the  path  was  an 
attempt  at  tightening  one  or  both  of  these.  In  a  1968 
paper,  Saltzburg  [4]  has  presented  a  tightened  upper  bound 
on  the  probability  of  error  for  a  linear  system  utilizing 
multilevel  pulse  amplitude  modulation  (M-ary  PAM)  having 
ISI  and  additive  Gaussian  noise.  Using  what  is  known  as  the 
Chernoff  bound,  in  which 

Pr[z  >  x]  <_  exp(-Ax)  E[exp(Az)]  X>0,  (39) 

where  E[y]  is  the  expectation  of  y,  a  bound  is  found  that  is 
the  equivalent  to  the  exponential  part  of  a  normal  distribution 
where  the  larger  ISI  components  act  to  reduce  the  signal  level 
and  the  smaller  components  add  to  the  noise  level,  with 
the  basic  system  of  Figure  4  being  considered,  the  PAM  signal 
input  to  the  demodulator  can  be  thought  of  as  the  amplitude 
of  the  inphase  component  of  an  M-ary  PSK  signal.  In  other 
words 


r<mTs)  =  a^lOI  +  k/0  am_kp(kTs)  +  nc(mTs) 

=  V’(0)  +  k/O  2k  +  nc(mTs) 


=  amP(0)  +  Z 


where  p(t)  represents  the  convolution  of  the  transmitted 
pulse  with  the  bandlimiting  channel  filter,  n  is  the  thermal 

v 

noise,  and  the  summation  in  (41)  represents,  of  course,  the 
ISI .  If  the  sequence  {a^}  has  the  M  equally  probable  values 


■j  -  2i  -  M  -  1 


j  —  1,2, 


then  the  overall  symbol  error  probability  is  given  by 

PJe)  =  Pr[z  >  p-] 


with  p£  -  |p(kTs)|.  The  Chernoff  bound  on  the  probability 
of  error  is  presented  as 


PM(e)  exp  " 


tPo-^-^keK  pk]2 

[2(an)2+j(M+1)  (M-D  kJK(p£)2] 


2  Z 

where  on  is  the  variance  of  nQ,  keK  p£  <  is  a  working 


wuclc  Uyj  -L.O  uuc  vaiiaucc  ±  llQr  keK  ®  wwj.rtAU'j 

constraint,  and  keK  represents  those  interferers  that  are 
assumed  to  subtract  from  the  signal  while  k^K  represents 
those  that  add  to  the  noise.  After  presenting  this,  Saltzburg 
gives  a  method  for  optimizing  the  set  K  for  the  minimization 
of  the  bound.  Finally,  (44)  is  applied  to  the  case  of  ideal 
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bandlimited  signaling  with  bandwidths  at  the  Nyquist  level 
and  below  by  way  of  example. 

Following  the  concepts  mentioned  above,  Lugannani  [5] 
presents  a  variation  on  the  development  of  the  error 
probability  expression.  He  recognizes  that  a  weakness  of  the 
Chernoff  bound  is  that  it  is  progressively  more  inaccurate 
for  small  boundry  values,  i.e.,  for  small  values  of  x  in 
?r[z>_x].  Consequently,  the  error  expression  is  divided  into 
two  terms  only  one  of  which  has  a  Chernoff  bound  applied, 
and  that  term  has  a  factor  that  is  small  when  x  is  small  thus 
reducing  the  uncertainty  introduced  in  the  bounding  process. 
Additionally,  the  bound  used  is  computationally  simpler  than 
that  used  in  (44)  yet  is  constrained  to  always  being  no  worse 
than  the  worst  case  bound.  The  system  model  used  by  Lugannani 
is  the  same  as  that  used  by  Saltzburg  with  the  exception 
that  M  *  2  in  order  to  simplify  the  analysis.  Due  to  the 
separation  of  the  effects  of  the  ISI  and  the  noise,  the 
demodulator  input  at  time  t  =  mTs  is  written 

r(mTs>  =  amp(0)  +  k£„  am_k  p(kTs)  +  nc(mTs) 


=  a  p{0)  +  z  +  n  (t) 
m  me 


and  the  expression  for  the  probability  of  error  evolves  into 


P,  (e)  *  h  Prl  |n  (t)  |  >  p'] 


+  *5  fQ  Iq(pQ-s)  -  q(Po+s)  ] 


Pr [ I z  >  s]ds 

i  m 


where  q(*)  is  the  pdf  of  n  .  If  it  is  assumed  that  B(0) 

^  2 

represents  the  worst-case  ISI  and  o  is  the  variance  of  the 
ISI  term,  then  the  bound  presented  in  this  paper  can  be 
derived  from 


PrUzm!  >  s]  <  min  [2exp (-s2/so,2)  ,  x(s)] 

2 


with  x(s)  defined  by 


(1  0  <.  s  <  B  (0) 

x  ( s )  *  j 

'  0  B  (0)  <  s 


Given  the  following  definitions, 


and 


n  =  min  [ X ,  In  4 ] ,  N  =  max  [X ,  In  4 ] , 
u  =  (1  -  y2)'H 


P2(e)  1  nYU$(s)ds  +  *  /"+nvV}(s)ds 

+  v  exp[-J5n2(l-Y2)]  fl*  \  |J^Ms)dS 
-  y  exp[-4n2  (1-y2)  1  f^u  +  nv5(s)ds 


where  4>(s)  is  the  normalized  Gaussian  density  function. 
Following  the  derivation  of  the  bound  in  (53)  ,  three 
different  pulses:  An  ideal  bandlimited  pulse,  a  fourth- 
order  Chebyshev,  and  a  Gaussian  are  used  as  examples  to 
compare  the  results  of  (53)  to  those  using  worst-case  bounds 
and  to  the  truncated  pulse  train  method  similar  to  that  of 
[3].  In  comparison  with  [4],  (53)  is  seen  to  be  tighter 
in  all  cases,  with  the  possible  exception  of  those  where 
the  Nq's  are  very  small  and  K  is  optimized.  In  addition, 

(53)  is  considerably  less  arduous  to  compute. 

The  next  logical  step  in  developing  an  upper  bound  on 
the  PM(e)  due  to  ISI  and  additive  noise  was  an  extension  of 
the  bounding  techniques  of  [4]  and  [5]  to  M-ary  PSK  modulation 
schemes.  Dealing  only  with  coherent  PSK,  Prabhu  [6]  applied 
the  methods  and  equations  of  [4]  to  binary  and  quarternary 
PSK  signals.  In  the  former  case,  he  notes  that  the  modulation 
is,  effectively,  only  applied  to  one  of  the  two  orthogonal 
carriers  and  thus,  the  situation  is  the  same  as  that  in  a 
binary  PAM  system.  The  bound  arrived  at  for  BPSK  is  therefore 
the  equivalent  of  that  given  in  equation  (44) .  Finding  a 
bound  for  the  QPSK  case  is  a  more  complicated,  albeit  very 
slightly,  task  than  that  for  BPSK.  With  the  (normalized) 
desired  signals  on  the  inphase  and  quadrature  carriers 
defined  as  aQ  and  bQ  respectively,  and  with  the  (normalized) 
combined  ISI  and  noise  for  the  two  denoted  by  za  and  zb 
respectively,  the  error  probability  for  a  given  phase 
$  =  $k  can  be  written  as 

P4(e)  =  hPr[a0+za  <  cos  <(>k]  +  *sPr[b0+zb  <  sin  $k] 

-  Pr[a  +z,  <  cos  <)>,  and  b  +z.  <  sin  ] 

o  a  k  ob  *  (53) 

=  Pr[(ao+za,  bQ+zb)eR1UR2] 

-  Pr[(aQ+za,  b0+zb)eR1DR2] 


....  A.-,.  J.-E. 
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where  and  R2  are  regions  of  probability  space  as  seen  from 
Figure  5.  Due  to  the  intuitively  small  probability  of 
R-jHr^  Prabhu  applies  the  "union  bound" 

P4(e)  <  Pr[(aQ+za,  b0+2b)£R1]  +  Prt(ao+za,  b0+zb)sR2]  (54) 

in  a  move  that  can  be  seen  to  reduce  (53)  to  a  sum  of  two 
BPSK  bounds — one  for  each  carrier.  Following  the  development 
of  the  bounding  equations,  several  graphs  are  presented 
showing  the  performance  of  various  pulses  as  a  function  of 
the  number  of  poles  in  an  RF  filter  and  the  filter  bandwidth. 

An  extension  of  Lugannani's  error  bound  to  M-ary  PSK 
systems  is  made  in  [7].  Following  the  same  developmental 
process  as  in  [5],  Koubanitsas  splits  the  error  expression 
in  two  parts  and  applies  the  same  Chernoff  bound  to  the  term 
representing  the  effects  of  ISI.  Since,  for  general  M-ary 
systems,  one  always  encounters  the  form  of  the  error  probability 
expression  as  given  in  the  paragraph  above  for  the  QPSK  case, 
the  union  bound  has  been  applied  for  simplification  of  the 
computations.  The  final  result  is  the  sum  of  several  error 
and  co-error  functions.  There  are  several  graphs  presented 
comparing  the  bound  that  was  derived  with  various  other 
bounds  including  that  of  [5]. 

Using  the  Chernoff  bounding  theory,  however,  is  not  the 
only  approach  that  has  been  taken  towards  finding  the  effects 
of  ISI  on  the  BER  of  a  linear  channel.  In  a  1974  paper, 

P.  J.  McLane  [8]  developed  a  lower  bounds  on  PM(e)  that,  when  used 
in  conjunction  with  any  of  several  previously  developed 
upper  bounds,  specifies  the  degradation  due  to  ISI  and  noise 
to  a  finer  level  of  precision  than  that  gained  in  using  an 
upper  bound  alone.  Considering  the  same  system  as  [4]  and  [5] 
with  the  ISI  assumed  to  remain  effectively  constant  beyond 
some  finite  number  of  pulses,  N,  and  to  have  a  symmetric  pdf. 


then  the  input  to  the  demodulator  at  t  =  mTg  is  written 

N/2 

r (mTc)  =  amp ( 0 )  +  D  avp(kTc)  +  nc(mTs) 


k=N/2 


=  amp(0)  +  +  nc(mTs) 


(41c) 


N  i 

with  the  IS I  from  all  pulses  beyond  j^|  assumed  to  be  0. 

The  denotes  summation  with  the  exclusion  of  the  k=0  term. 

XI 

If  it  is  assumed  that 


Gn(x)  =  Pr(nc  >_  x) 


=  1  -  Fn(x) 


where  F  (x)  is  the  distribution  function  of  n  .  then  the 
n  c 

expression  for  the  error  probability  may  be  given  in  the  form 
P2(e)  =  <i  (ez  [Gn(p-  +  |zm|)]  +  Ez  [Gn(p'  -  Uml)l)  (56) 


withEy[g(y)]  denoting  the  expectation  of  g(y)  with  respect 
to  Y.  At  this  point,  a  final  condition  is  placed  on  the  signal 
for  xeX,  Gn(x)  is  assumed  to  be  convex.  An  application  of 
Jensen's  inequality  is,  therefore,  allowed  and  yields 

P2(e)  >  MG„tPo  +  E(!zm|l!  +  Gn[Pp  "  E(lzml)!  )  (57) 


for  all  p  ±  E(|z  I)  e X.  Since,  in  general,  E ( ] z  ! )  is  almost 
o  ‘  m  m 

impossible  to  compute,  upper  and  lower  bounds  for  this  term  are 
needed.  The  upper  bound  is  shown  to  be  a-  whereas  two  lower 


bounds,  b^,  are  given  as 


where  the  use  of  yields  what  is  known  as  the  "fourth 
moment  lower  bound"  and  bz  the  "ratio  lower  bound."  McLane 
proceeds  to  present  several  graphs  showing  his  bounds  and 
those  of  some  others  for  various  signal  pulses  and  values 
of  N. 

A  developmental  branch  of  the  analytical  methods  that 
arose  concurrently  with  the  bounding  approach  uses  the  Gram- 
Charlier  expansion  of  the  ISI  (or  ISI  +  noise)  pdf  to  allow 
direct  evaluation  of  P(e)  to  just  about  any  desired  level  of 
accuracy.  In  a  1971  paper,  Shimbo  and  Celebiler  [9]  first 
introduced  this  method,  applying  it  to  binary  PAM  signals  in 
the  system  of  Figure  4.  In  a  parallel  to  the  work  of  Saltzberg 
[4],  the  particular  approach  used  involves  the  expansion  of 
the  pdf  of  the  noise  +  ISI.  The  demodulator  input,  therefore, 
has  the  form  of  (41a)  for  t=mT  and  the  probability  of  error 

w 

can  be  written  as 


P,  (e)  »  [1  -  Pr  ( |  z  |  £  p')  ] 


(60) 


=  Ml  -  Qe). 


Qe,  as  given  above,  can  be  expressed  in  terms  of  the  distribution 
function  of  Z ,•  f2(z),  which  in  turn  can  be  expressed  as  a 
Fourier  transform  of  its  characteristic  function.  Since  the 


characteristic  function  of  the  ISI  and  the  AWG  noise,  Q 
is  given  as 


e' 


(61) 


Qe  =  FZ(po}  '  FZ("po) 


l  r=eSin(p'u) 
7  Jx  u 


2  2 

(H  cos  p/u)  exp  (=*50  u  )  du 
^  jc  n 


! 


V 


where  the  infinite  product  of  cosines  is  the  characteristic 

function  of  the  ISI  and  the  exponential  factor  is  that  of  the 
2  . 

noise.  0n  is  the  noise  power.  If  the  following  definitions 
are  allowed. 


2 

a 


+ 


s=Pq/o , 


Vpk/fl' 


(62) 


the  Qe  can  be  rewritten  as 


Q 


e 


1  f  <*>  sin  (su) 

it  J  u 

QO 


exp  (hau 


2 

(II  cos  a.  u)  exp(-J$u  )du. 

1  K 

(63) 


The  Gram-Charlier  power  series  expansion 


2  2k 

exp(*sau  )  (n  cos  a.  u)  =  1  +  Z  b,vu 
1  k  l 


(64) 


can  now  be  applied  and  the  resulting  expression  for  the  error 
probability  is 


P2(e)  =  erfc(s/2)  +  I  (-1)  b2kG2k-l^s* 


(65) 


when  Gj^  ( s )  is  the  k-th  order  Hermite  function 


V>  -  ^4  4  <«*p  -*2> 


O.  ■OTV  ^ 


GJc+ 1  ( s )  =  sGk(s)  -  kGk_1(s) 


(66) 


The  coefficients  of  the  series  expansion, 
be  expressed  as  the  recurrence  relation 


b,.  ,  are  found  to 


1  k 

b2k  “  2k  E.  b2k-2nd2n-l 
n=l 


(67) 


where  bQ  =  1  and  d2n-1  is  the  product  of  the  (2n-l)-th 

coefficient  of  the  power  series  expansion  of  tan  u  by  the 

sum  E^»ak^nj  =  0*  In  equation  (65),  ?2(e)  consists  of 

two  terms,  the  first  of  which  can  be  thought  of  as  the  error 

2 

due  to  additive  Gaussian  noise  of  power  a  and  the  second  as 

that  due  to  the  infinite  IS I.  If  the  intersymbol  interference 

2  2 
power,  is  smaller  than  the  noise  power,  an  ,  then 

(65)  will  converge  rather  quickly.  If,  however,  the  ISI  is 

large  compared  to  the  noise  power,  convergence  will  be  slow, 

and  thus,  another  expression  for  p2 (e)  is  given  that  causes 

the  rate  of  convergence  to  increase: 

p2(e)  =  j-  (  E  lh  erfc(si/2)  +  E(-l)  kb2kG2k-l  (sj>)  ])  (68) 

'i-1  2  ' 


with 


si  =  Pi/t',n2+j=r(Pj)2l'!  1  “l'  1  “2' 

V  Pk/[VS.r<Pj>2l!i 


a 


r- 


1 


(69) 
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This  means  that  the  error  calculation  becomes  the  enumeration 
of  the  ISI  for  r-1  terms.  A  method  of  calculating  the 
resulting  truncation  error  is  presented  and  (68)  is  applied 
to  various  pulse  shapes  as  a  function  of  L  and  compared  to  some 
previously  developed  P2(e)  bounds. 

Once  again,  the  next  logical  step  in  the  development  of 
the  series  approach  was  an  extension  of  the  work  of  [9]  to 
M-ary  PSK  (CPSK)  modulation  schemes,  and  once  again  Prabhu 
was  the  one  to  perform  such  an  extension.  In  [10]  such  an 
extension  is  made.  After  passing  over  the  binary  case  with 
just  a  few  brief  comments,  the  case  of  quarternary  PSK  is 
chosen  as  the  primary  vehicle  for  development  of  the  bounds. 
Utilizing  the  same  demodulator  input  signal  as  in  [6] ,  an 
expression  for  the  error  probability  given  that  ^  J  at  t  =  t 
is 


P(e  |*k«J)  =  1  -  Pr[a(to)+za(tQ)  ^0  and  b  (tQ)  +  zb  ( tQ)  >_0]  (70) 


If  it  is  now  assumed  that  the  major  interference  terms 
correspond  to  -N]_  <  M  <  N2 ,  N  =  Nj.  +  N2 ,  then  the  conditional 
probability  can  be  written 


P  (e  I 


>k-T)  = 


1  -  Pr[yn+yR  >  y 


r0C 


and 


ZN+zR  -  "  Yos^ 


(71) 


=  W  +  meA^m 


where 


zXT  =  s  ( t  )  +  Z.t_ 
N  o  meA  in 


(73) 


ZR  5  m l 


ATnT  ECzr^  "  0 


A  =  a  (y  )  -  E [a  (t  )  ] 
m  m  ■*  o  mo 


5  b  (  t  )  -  E  [  b  ( t  )  ] 
m  m  o  mo 


(74) 


and 


A  =  a (t  )  +  m5n  E[am(t  )] 
oc  o  m?*o  m  o 


xoc  =  b (t  )  +  m5n  E[b  (t  ) ] 
os  o  m^o  m  o 


(75) 


From  this  point,  it  is  shown  that  the  bounds  for  the  probability 
of  error  can  be  written  as 

F„  (-y„„-M  +  P,  (-y«.*A)  -  P„  „  A,-  og+A) 


yN  oc 


ZN  03 


yNZN  OC 


-Pr[|yR|>M  -  Pr  [  j  zRj  >  a]  P  (e !  4>lc=J)  -FyN(“Yoc+A) 


+  fzn(‘tos+1)  -  ’’'os-41 


(76) 


+  Pr[ |yR! >a]  +  Pr [ | zR| >a] 

where  F^(x)  =  Pr[y  <_x]  and  A  is  any  real  number.  Following 
this,  Prabhu  uses  the  Gram-Charlier  expansion  to  compute  the 
F(*)'s  and  proceeds  to  present  a  method  of  bounding  the 
error  whereby 


fl(meA)  -  5(A)  <  P(e)  <  Q(meA)  +  5(A)  (77) 

with  fl(*)  the  series  expression  of  the  probability  of  error 
due  to  the  IS I  terms  where  me  A  and  5(A)  is  the  truncation 
error  for  those  terms  me A.  Bounds  for  5(A),  and  hence  for 
P(e),  are  shown  to  be  computable  using  the  Chernoff  technique 
and  the  optimization  of  A  is  performed  the  same  way  as  in 
[4].  An  extension  from  QPSK  to  general  M-ary  PSK  is  discussed 
and  several  graphs  and  diagrams  are  given  showing  the 
performance  of  this  approach. 


3.1.2 


Approaches  to  Nonlinear  Systems  with 
Additive  Noise  Only 


As  the  number  of  satellites  and  satellite  users  grew, 
it  rapidly  became  apparent  that  the  power  and  signal-bandwidth 
limitations  inherent  in  FDMA  were  bound  to  doom  the  technique. 
TDMA  was  obviously  the  solution  to  the  problem,  yet  the 
operation  of  the  satellite  transponder  to  achieve  peak  power 
efficiency  created  its  own  problem:  a  nonlinear  TWTA  transfer 
characteristic.  It  became,  therefore,  necessary  to  determine 
the  effect  of  this  nonlinearity  on  the  system  performance. 

The  next  level  of  complexity  in  the  approaches  to  the  performance 
prediction  problem — above  and  beyond  those  dealing  with  the 
effects  of  ISI  and  additive  noise  in  a  linear  system — are 
those  attempts  at  the  determination  of  the  effects  of  additive 
noise  in  a  nonlinear  system.  Obviously,  some  of  the  noise 
must  be  assumed  to  appear  at  the  input  to  the  TWTA  (i.e.,  the 
uplink  noise)  and,  indeed,  all  of  the  work  summarized  herein 
makes  that  assumption.  Descriptions  of  the  various  approaches 
to  this  problem  follow.  It  should  be  noted  that  various  and 
sundry  reasons  are  given  for  ignoring  the  existence  of  ISI  and 
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that  these  reasons  will  be  mentioned  within  the  summary  of 
each  approach. 

In  Figure  6,  a  model  of  a  nonlinear  system  with  additive 
noise  is  shown.  As  could  be  surmised,  it  is  almost  identical 
to  that  of  Figure  3  with  the  bandpass  nonlinearity  representing 
the  whole  transponder  and  without  the  bandlimiting  character 
of  the  transmit  and  receive  filters.  Given  that  this  system 
has  the  same  modulator  output,  D  ,  as  does  the  system  of 
Figure  3,  then  the  input  to  the  bandpass  nonlinearity  is 


where  n  is  the  same  as  that  of  (23)  and  p  is  the  same  as  that 
in  equation  (17).  For  a  particular  time,  kTs  <_  t  <  (k+l)Ts, 
it  can  be  seen  that  there  is  no  ISI  present.  Passing  through 
the  bandpass  nonlinearity,  once  again  described  by  f(*)  and 
g('),  the  signal  at  the  input  to  the  receive  filter  can  be 
given  as  in  (28)  exclusive  of  the  effect  of  the  ISI  introduced 
by  the  transmit  filter  of  Figure  3.  Since  there  is  also  no  ISI 
as  created  by  the  receive  filter,  the  demodulator  input  for 
the  time  kTs  <  t  <  (k+l)Ts  can  be  expressed  by 


FJgure  6  -  Bandpass  Nonlinearity  with  Noise 


Pk  cos ( t  +  ek) 


(81) 


Y_  = 


where  ok  is  as  written  in  (33)  (excluding  the  ISI  that 
affects  Wk  as  was  previously  noted) . 

Although  some  earlier  work  has  been  done  on  the  system 
of  Figure  6,  one  of  the  first  approaches  to  the  problem  was 
presented  in  a  1972  paper  by  Davisson  and  Milstein  [14], 

This  approach  comprised  two  parts,  the  first  concerned  the 
case  of  a  one-link  system  where  the  downlink  of  Figure  6  is 
excluded  and  the  second  being  concerned  with  a  two- link 
system  similar  to  that  in  Figure  6  with  a  second  bandpass 
nonlinearity  at  the  receiver  input.  All  of  the  nonlinearities 
are  assumed  to  be  identical  limiters  having  the  transfer 
characteristic  shown  in  Figure  7.  Additionally,  both  the 
one  and  two-link  systems  are  assumed  to  have  a  demodulator 
output  that  is  the  result  of  shifting  the  signal  down  to 
baseband,  integrating  for  T  seconds,  and  sampling. 

In  the  first  part  of  (14],  the  modulation  scheme  given 
is  BPSK  and  the  input  to  the  bandpass  limiter,  therefore,  is 


Xp  =  (Rk  cos  *k*  cos(u)Qt  +  6k) 
=  Vki  ( t)  cos  (u)Qt  +  6k) 


(82) 


where  all  of  the  variables  are  as  previously  defined  and 
$k  *  0  or  it.  Since  the  limiter  only  has  an  AM- AM  character¬ 
istic  which  we  can  describe  by  f(*),  then  the  output  of  the 
limiter  is 


Zp  =  f(Vki}  cosUot  +  skJ 


=  vko(t)  cosUot  +  Sk) 


(83) 


50 


*  -*  ’•*  m 


and  the  sampler  output,  d(t),  can  be  expressed,  for 

t=(k+l)T  ,  as 
s 


(k+1)  T 

d[(k+l)Ts]  =  h  J  \0  cos  3k  dt*  (84) 

kTs 


It  is  noted  that  (84)  is  a  stochastic  integral  that,  in 
general,  has  an  unknown  pdf.  To  make  the  theoretical  analysis 
tractable,  the  baseband  detector  is  henceforth  assumed  to  be 
of  the  sample-and-sum  type  instead  of  the  integrate-and-dump 
as  was  described  previously.  If,  further,  k=0,  this 
assumption  yields 


d(Tg) 


(85) 


With  the  noise  power  spectral  density  assumed  to  be  flat  over 
a  bandwidth  equal  to  2  ( 2 tr f  ),  N  is  chosen  to  be  2f  T  in  a 

C  w  O 

move  to  insure  the  independence  of  the  samples  in  (85) .  At 
this  point,  the  term  to  the  right  of  the  l  in  (85)  is  defined 
as 


a 


i 


(86) 


and  the  error  probability  is 

r2f  _T 

I  / 

P2(e)  = 


/2fcTs 

H  Prl  £  < 

'  i  =  1  1 


0  4> 


/2fcTs 

+  Pr^  £  aj_  >  0|  p 
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(87) 

] 


■  ^  ji*'  *"*■ 


from  which  it  can  be  seen  that  the  calculation  of  P2(e) 
requires  that  the  probability  density  of  a.  be  known.  This 
pdf  is  derived  and  is  given  by 


q  (a)  = 


where 


and 


The  derivation  of  (88)  should  be  noted  as  being  somewhat 
complicated  and  evaluation  of  P2 (e)  is  performed  by  using 
a  computer  to  perform  the  necessary  2fcT  numerical  convolutions 
of  (88) .  A  graph  is  presented  showing  the  relative  performance 
of  the  system  as  calculated  from  (87)  as  a  function  of  the 
e/pQ  ratio  and  in  comparison  to  the  "continuous-time  linear" 
system  performance.  It  is  seen  that  the  closer  the  limiter 
comes  to  being  an  ideal  hard  limiter  (e=0) ,  the  more  the  one 
link  system  performance  is  degraded — up  to  2  dB  for  a  BER  of 
1Q-^  in  the  ideal  case. 

For  the  two-link  system  previously  mentioned,  Davisson 
and  Milstein  are  effectively  just  repeating  the  work  of  the 
first  part  of  their  paper.  Although  they  discuss  a  variation 


o  >  1 


2  e'A(e'ct-X)  [erf  (1-a2)  ') -h] 


+  (2/ir)  ^expf-^A2  (1-a2)  /  (l-ct2)*5]  [  ( Xa-e  ") 


(88) 


L 


+  Xaerf (Xa-e  ') ] 


a!  <  1 


X  = 


*k 


°n2  ' 


(89) 


n 


A  (x)  =  (2ir)  55  exp(-x2/2). 


(90) 


of  the  two-link  system  that  was  studied  by  I.  Jacobs  in  a 
1965  paper,  it  is  a  very  limited  discussion  even  though  there 
is  an  extensive  expansion  and  solution  of  the  error  equation 
given  by  Jacobs  presented  in  an  appendix.  The  sum  total  of 
this  study  of  a  two-link  system  with  two  bandpass  limiters 
was  that  graphs  were  generated  comparing  the  performance  of 
the  authors'  system,  the  system  studied  by  Jacobs,  and  a 
similar  linear  system  for  three  different  values  of  N  =  2fcTs, 
the  number  of  samples  taken.  Resulting  from  these  comparisons 
is  the  conclusion  that  multiple  limiters,  with  and  without 
Jacobs'  clipper,  enhance  a  system's  performance,  particularly 
for  small  values  of  N,  over  that  expected  for  linear  systems 
with  the  same  number  of  noise  sources  (>_2)  .  This  is  in  contrast 
to  Jacobs  who  predicts  enhancement  only  for  small  BTg  products. 

Following  the  developmental  process  initiated  in  [14] , 
one  could  presume  that  the  next  logical  step  would  have  been 
the  extension  of  the  nonlinearity's  characteristic  from  that 
of  a  limiter  to  a  more  general  case.  Additionally,  one  might 
also  expect  an  extension  of  the  modulation  format  to  a 
generalized  M-ary  PSK  one.  This  was  indeed  the  situation  in 
a  1973  paper  by  Lyons  [15].  Using  analog  PM  notation  in 
analyzing  the  effect  the  nonlinear  transponder  has  on  a  signal 
(an  acceptable  technique  that  is  notationally  convenient  and 
can  be  generalized  with  ease  to  the  PSK  case  since  there  is 
no  bandlimiting) ,  the  transponder  output  is  the  same  as  Z^ 
of  Figure  6.  Lyons,  however,  separates  Z  into  the  desired 
and  undesired  portions  of  the  signal.  Recalling  that  (using 
analog  PM  notation) 

Z  =  f(R)  cos[oj  t  +  <j>  +  y  ( t )  -  g(R)].  (91) 

P  o 

where  R  and  i>  are  the  analog  representations  of  the  previously 
defined  ((78) -(80)]  variables  RK  and  and  y(t)  is  the 
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phase  error  due  to  the  uplink  noise.  The  desired  portion 
of  the  signals,  Z'  (t) ,  is  found  by  taking  the  expectation 

Jr 

of  Zp  conditioned  on  the  known  value  of  u>Qt  +  <j>.  This  yields 


wi  th 


A  cos (uQt  +  $  -  n) 


(92) 


A  =  E ( f  ( R)  cos  [ y  -  g(R)  +  n ] )  (93) 


and 


n 


arctan 


/ E  (f  (R)  [sin  y  ~  g  (R)  ])\ 
\E  (f  (R)  [cos  y  -  g  (R)  ])/ 


(94) 


Obviously,  therefore,  the  unwanted  portion  of  the  signal,  the 
"pseudo-noise,"  n"(t),  may  be  written  as 


n'  (t) 


Z 


P 


n'  (t)  cos(u  t+<j>-n) 
c  o 


n'  (t)  sin  (<*>  t+$-n) 
s  o 

(95) 


nc  (t) 

=  f  (R) 

COS  [ y 

-  g  (R)  +  n]  -  A 

(96) 

ns  (t) 

=  f  (R) 

sin  [  y 

-  g (R)  +  n] . 

(97) 

The  downlink  noise,  v,  now  combines  additively  with  the 
"pseudo-noise"  and  the  path  that  needs  to  be  taken  becomes 
pretty  clear  albeit  exceedingly  complicated.  To  reduce  the 
complexity,  Lyons  makes  the  assumptions  that,  1)  a  number  of 
samples  of  the  received  phase,  N,  are  taken  over  one  symbol 
duration  and,  as  in  [14],  to  insure  their  independence 
N  =  2n+1^2fcTs,  and  2)  the  detector  operates  by  majority 
decision  logic.  Taken  together,  these  two  assumptions  make 
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the  error  probability  the  probability  that  n  +  1  samples 
are  incorrect;  this  is  given  by 


ve>  - 


f2,ri')[i-i 


j0rnu-p 


j^(e)  ]  [P^(e)  ] 


2n+l-k 


(98) 


where  PM(e)  is  the  probability  of  error  for  a  single  sample 


Pr[  |  6  | 

>  ff/M] 

(99) 

f  n'  +  v 

arctan 

U  +  V  A  ] 

(100) 

To  solve  (99) ,  it  is  necessary  to  know  the  joint  pdf  of 

n",  n',  v' ,  and  v  .  Since  n  and  v  are  assumed  to  be 
esc  s 

statistically  independent  with  known  pdf's,  then  the  solution 

can  be  achieved  just  by  calculating  q(n',  n') .  A  general 

method  of  doing  this  is  given.  An  example  and  related  graphs 

are  presented  for  the  case  in  which  the  bandpass  nonlinearity 

is  a  limiter  (similar  to  the  limiter  of  Davisson  and  Milstein's 

paper) .  The  graphs  are  given  as  functions  of  the  signal-to- 

noise  ratio  (SNR)  and  the  B  T  product.  As  in  [14],  the 

n  s 

limiter  is  seen  to  degrade  performance  in  the  absence  of  down¬ 
link  noise  and  enhance  performance  in  its  presence  when  (SNR)  = 
(SNR)down  for  the  BPSK  case — although  not  for  QPSK  or  3-PSK. 
Overall,  Lyons  notes  that  the  system  degradation  due  to  a 
hardlimiting  nonlinearity  is  an  increasing  function  of  M,  the 
percentage  of  the  total  noise  that  is  found  in  the  uplink,  and 
the  BnTg  product,  except  for  BPSK  when  the  downlink  noise 
dominates  in  which  case  the  performance  is  enhanced. 

Extending  the  work  done  in  [14]  and  [15],  Jain  and 
Blachman  [16]  look  at  the  system  of  Figure  6  with  a  BPSK 
modulation  format  and  an  ideal  hard  limiter  as  the  system 
nonlinearity.  The  primary  intent  of  their  study  was  to  justify 
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and  demonstrate  the  desirability  of  limiting  via  an  accurate 
method  of  calculating  the  error  performance.  Making  the  same 
assumptions  vis'  selection  of  the  number  of  samples  taken 
per  bit  as  we  have  discussed  previously,  the  error  probability 
expression  becomes  that  presented  earlier  by  Lyons  (98) . 

Since  the  possible  values  of  the  modulation  phase  are,  as 
usual,  considered  to  be  equally  probable,  the  calculation  of 
Pm(E)  devolves  into  the  calculation  of  the  error  probability 
PM(e)  for  any  one  sample  of  the  received  signal.  Due  to  the 
ideal  nature  of  the  limiter,  where  f(*)  =  1  and  g(*)  =  0, 
the  transponder  output  is  written  (in  analog  PM  notation)  as 


Z  =  cos  (u>  t  +  i) 
P  o 


8 


arctan 


arctan 


r  posin  » +  ns1 

L PQcos  ♦  +  nc J 

r  ns  l 

I_p  cos  $  +  n  J 


(101) 


(102) 


and  the  receiver  input  as 

Yp  =  a  cos(u)Qt+8)  +  vc  cos  (i>Qt  -  vg  sin  u>Qt  (103) 

where  a  is  determined  by  downlink  losses,  etc.,  and  heretofor 
has  been  assumed  to  be  unity.  With  a  coherent  demodulator 
and  remembering  that  sin(rnr)  =  0  for  all  n,  the  sampled  base¬ 
band  representation  of  Y  becomes 

Jr 


r (t  )  =  a  cos  6 ( t  )  +  v  ( t  )  . 


(104) 


In  order  to  find  P2  (e)  it  is  now  seen  to  be  necessary  to 
determine  the  pdf,  q(r),  of  the  sampled  output.  From  the 
known  pdf  of  vc  and  (104)  this  is  obtained  by  convolving  the 
pdf  of  the  phase  8  with  that  of  the  noise  yielding 

q (r)  =  /027r  exp  ^(r-^CO|  B)^  q(6)  dS. 

v  (105) 


If  it  is  assumed  that  the  transmitted  bit  is  a  1  (4>  =  0), 
P2 (e)  can  now  be  written  as 


P2(e) 


=  f°  q  (r)  dr 

—  CO 


2*  r  2  ro  /  -(r-acos  s)2\  1 

=  SQ  <3  (8)  [(2avi 2)  *  fm  2q2  )  drj 


dB 


=  H  -  h  f^erf  (a[2  ( ov)  Z]~ii  cos  b)  q(3)  dB 


2,-Js 


(106) 


from  which  it  is  easy  to  see  that  P2(e)  is  unaffected  by  the 
limiter  when  all  of  the  noise  is  on  only  one  of  the  two  links. 

It  should  be  noted  that  this  conclusion  contradicts  [14]  for 
the  case  of  no  downlink  noise.  The  error  function  is  now 
expressed  in  the  form  of  a  Fourier  series  containing  an 
integral,  the  solution  of  which  is  given  in  terms  of  a  confluent 
hype rgeome trie  series.  Consequently,  the  probability  of  error 
is  expressed  by 


P2(e)  = 


1  -  i  E 
A  tt  n=( 

(n 


lFl\n+H,2n+2, 


<-l)nr(n+J5)  r.,o_  2  x  — *5 1  2n+l 

(2'n+l)  ! -  la(2av  }  J 

- ^o-V2lT  cos  [  (2n+l)  3  ]  q  (  3) 

2°  2J  ° 

i  V  ' 


(107) 
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where  ^F^  ( • )  is  the  confluent  hypergeometric  series  and 
r  ( • )  is  the  gamma  function  as  usual.  Following  the  derivation 
of  (107) ,  three  different  yet  equivalent  expressions  for 
P2  (e)  are  presented.  In  the  first,  an  integral  representation 
for  q(6)  (for  a  specific  value  of  <j>)  and  a  convergent  infinite 
series  solution  is  developed.  The  second  expression,  again 
a  series  solution,  is  conjured  up  by  using  a  relationship 
whereby  the  confluent  hypergeometric  series  is  a  function 
of  two  sequential  modified  Bessel  functions.  Lastly,  an 
equivalent  to  (107)  is  developed  through  the  use  of  a  closed 
form  representation  of  q(S).  This  leads  to  an  evaluation  of 
(106)  that  gives  P2  (e)  in  terms  of  Rice's  Ie  function 

Ie(k,x)  » /x  exp(-t)  I  (kt)  dt.  (108) 

o 

For  the  special  case  when  SNRu  =  SNR^  »  p^,  the  error 
probability  is  simply 

P2(e)  =  H  exp (-  P;L2).  (109) 

Jain  and  Blachman  proceed  to  give  numerical  results  comparing 
their  three  P2(e)  expression  vis'  ease  of  calculation  and 
accuracy  and  then  move  on  to  the  comparison  of  system  per¬ 
formance  with  and  without  the  limiter.  These  latter  comparisons 
indicate  that  the  limiter  produces  a  2-3  dB  performance 
enhancement  for  large  uplink  SNR  (>  10  dB)  and  a  somewhat 
smaller  magnitude  degradation  for  small  SNR  on  the  uplink. 

The  authors  conclude  their  work  by  deriving  an  optimum 
nonlinearity  to  enhance  system  performance. 

A  fairly  complete  picture  of  the  effects  a  limiter- 
type  nonlinearity  has  on  the  performance  of  wideband  satellite 
systems  with  AWGN  on  both  the  uplink  and  the  downlink  has 


been  presented  by  the  preceding  papers  [14] -[16].  Obviously, 
the  next  step  in  the  analysis  of  this  kind  of  wideband  system 
is  the  consideration  of  a  more  realistic  type  of  transponder 
nonlinearity.  This  step  was  taken  and  is  discussed  in  a  1976 
paper  by  P.  Hetrakul  and  D.  P.  Taylor  [17].  Using  a 
quadrature  model  of  the  TWTA,  the  authors  developed  a  model 
wherein  the  inphase  and  quadrature  AM-AM  characteristics, 
FjtRft)]  and  F^tlUt)]  respectively,  are  assumed  to  take  the 
form  (with  the  time  dependence  of  R(t)  taken  for  granted) 

Fr(R)  =  CjR  exp(-C2R2)  IQ  (C2R2)  (110) 

FQ(R)  =  SjR  exp(-S2R2)  I1  (S2R2)  (111) 

where  I  (•)  is  the  modified  Bessel  function  of  the  first 
n 

kind  of  order  n.  Using  any  TWTA  whose  transfer  characteristics 
are  known,  the  coefficients  of  these  equations  can  be  found 
through  the  use  of  the  fundamental  least-squares  curve¬ 
fitting  techniques.  For  illustrative  purposes  the  authors 
determined  the  coefficients  for  the  TWT  used  on  the  Intelsat 
IV  satellite;  equations  (110)  and  (111)  are  shown  to  yield 
very  close  fits  to  the  characteristics  of  the  TWTA  in  question. 
Representing  the  signal  in  the  manner  introduced  by  Lyons, 
Hetrakul  and  Taylor  present  an  expression  of  the  pseudo-noise 
pdf's  that  is  in  the  terms  of  a  Gram-Charlier  expansion. 

Using  the  notation  of  the  inphase  pseudo-noise  term, 


exp 


H 
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where  [ • ]  and  b2k 
quadrature  term,  of 


are  used  as  defined  in  [9]. 
course,  has  the  same  form. 


The 


(112) 


Following  the  determination  of  (112) ,  the  authors 
proceed  to  give  examples  of  BER  determination  for  BPSK  systems 
with  both  a  majority  logic  receiver  and  a  matched-filter 
receiver.  Performance  of  the  two  receiver  types  is  compared 
and  it  is  seen  that  the  quadrature  model  of  the  TWTA  that 
has  been  introduced  yields  results  that  agree  with  those  of 
Jain  and  Blachman  [16]  for  the  case  of  the  band-limiting 
nonlinearity.  Substantiating  evidence  for  the  results  is 
given  in  the  form  of  a  graph.  After  concluding  that  the 
performance  calculations  can  be  extended  to  the  M-ary  PSK 
system,  an  example  of  the  numerical  calculation  technique 
used  in  the  example  to  derive  the  moments  of  the  inphase 
noise  component  is  presented  along  with  proof  of  the  absolute 
convergence  of  the  derived  BER  expression. 

A  different  approach  to  the  consideration  of  a  more 
realistic  type  of  transponder  nonlinearity  can  be  found  in 
a  1977  paper  by  Forsey,  Gooding,  McLane,  and  Campbell  [18], 
Whereas  the  authors  of  [17]  developed  a  unique  quadrature 
model  of  the  TWTA,  Forsey,  et  al . ,  stick  with  the  previously 
introduced  method  of  using  the  generalized  methods  of 
representing  the  AM-AM  and  AM-PM  transfer  characteristics 
as  f(*)  and  g(*)  respectively.  Further,  the  previously 
discussed  paper  utilizes  the  "pseudo-noise"  representation 
of  the  transponder  output  that  was  derived  by  Lyons  [15]  while 
[18]  follows  the  analysis  techniques  of  Jain  and  Blachman 
in  [16].  In  fact,  the  work  presently  under  scrutiny  is  just 
a  generalization  of  the  analysis  presented  by  Jain  and 
Blachman  for  hard-limited,  BPSK  systems. 

The  authors  open  their  discussion  with  the  development 
of  the  signal  as  it  can  be  represented  at  the  input  to  the 
demodulator  with  the  end  product  being  the  same  as  that  given 
in  equation  (82) .  Following  a  rather  brief  description 
of  phase  compensation,  the  standard  error  equation  is  conjured 
with  the  form  for  the  case  of  the  transmitted  phase,  ^ ^=0 , 
being 
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tt/M 


qe|0k=o(9) 


(113) 


Solving  (113)  obviously  is  the  logical  progression  from  the 
determination  of  the  given  conditional  probability  density 
function  (hereafter  denoted  by  q  ( 6 ) )  and,  therefore,  the  meat 
of  the  work  done  in  this  paper  is  an  application  of  the 
analysis  of  [16]  to  that  determination. 

Forsey,  et  al.,  find  PM(e)  through  a  four-step  procedure 
that  begins  by  conditioning  on  R(t)  and  8(t)  and  drawing 
the  resultant  pdf 

q[e|R(t),  Bit)]  -  jV  ^  (HlPl(l.  n+1:  -V2) 

\  (114) 

cos[n(e  -  8(t)-g[R(t)  ])j  J 


.  2  _  fz[R(t)] 
Y2  "  ~T~2 


(115) 


from  elsewhere  in  the  body  of  current  knowledge.  r(*)  and 
(•;•;•)  are  the  gamma  function  and  the  confluent  h^per- 
geometric  function  respectively  and  have  been  defined 
previously.  Step  two  consists  of  solving  (113)  for 
q [ 0  j R(t) , 0 (t) ] — an  action  yielding  a  rather  long  and  complex 
result  that,  once  again,  has  been  taken  from  other  literature 
The  next  step  in  the  solution  process  is,  obviously,  the 
removal  of  the  conditioning  on  R(t)  and  6(t)  which  yields 

r°°  r  2* 

PM(e)  ■  /  /  PM[elR(t)'  S(t)]q[R(t),  6 (t) ]  dS(t)dR(t). 


Since  the  probability  density  above  is  just  the  joint  density 
of  the  amplitude  and  phase  of  the  signal  that  results  when 
narrowband  Gaussian  noise  is  added  to  a  sinusoid  of  known 
amplitude  and  phase — a  very  well-known  pdf — then  the  final 
step  of  the  solution  is  to  solve  (116).  The  resulting  error 
expression  is 


PM(e) 


=  M-l  +  (2exp  (~  Yn  )  5  r  (-1)  nr  (n/2+l)  1 

M  \  1ran2  n=1  Lnr(n+l)(/2av)nJ 

sin  ("-’ft— ’■)  5<nlj 


(117) 


(118) 


/  2 

5(n)  =  J  R(t)fn[R(t)  ]cos|ji(g[R(t)  ]  -g(A))J  expf-  R 

°  '  2an 

iFl(?'-n+1'-  «v2)  (-^)  dR(t)  (U9) 


where  In(*)  is  as  previously  defined.  It  is  stated  that  the 
infinite  series  in  (117)  can  be  shown  to  be  convergent.  A 
little  consideration  will  reveal  that  this  error  expression 
is,  indeed,  the  generalized  form  of  that  presented  in  [16] 
for  the  hard-limited  BPSK  case. 

Following  the  presentation  of  (117)  the  authors  proceed 
to  verify  their  work  through  the  use  of  an  example  utilizing 
a  TWTA  model  drawn  from  the  literature.  Although  no  numerical 
verification  was  presented,  it  is  shown  that  (117)  is  usable. 


After  a  discussion  of  the  effects  of  TWT  input  backoff 
on  system  performance,  the  work  is  concluded  with  the 
assertation  that  system  IS I  can  be  taken  into  account  using 
the  methods  already  described  by  Prabhu  [6],  [10]  in 
combination  with  basic  enumeration. 

The  work  presented  in  [18]  is  by  far  the  most  complete 
treatment  of  the  performance  prediction  problem  as  it 
pertains  to  wideband  nonlinear  satellite  communications 
systems  with  significant  uplink  noise.  Unfortunately,  when 
M  >  4  it  becomes  a  fairly  complex  task  to  evaluate  the  infinite 
series  that  appears  in  the  final  error  expression.  In  [19], 
authors  Mathews  and  Aghvami  have  recognized  this  deficiency 
and  present  an  approach  designed  to  skirt  the  complexity  of 
evaluation  while  retaining  the  accuracy  for  M-ary  systems 
where  M  >  4. 

Following  the  usual  signal  development  up  to  and  including 
the  input  to  the  demodulator,  the  authors  assume  that  the 
signal  is  coherently  demodulated  (with  phase  compensation  e) 
and  sampled  once  every  symbol  duration  yielding, 

I  =  f (R)  cos  [  B  -  g  (R) +e  ]  +  vc  (120) 

Q  =  f (R)  sin [ B  -  g (R) +e ]  +  vg  (121) 

as  the  inphase  and  quadrature  components  respectively. 

Using  very  generalized  reasoning,  the  basics  of  error 
probability  are  given  and  the  expressions  for  the  error 
probabilities  of  the  binary  and  quarternary  cases  conditioned 
on  R  and  8  are  developed.  After  stating  that  it  is  not 
possible  to  obtain  a  simple  analytical  expression  for 
P  (e|R,  B)  for  M  >  4,  the  original  contribution  of  [19]  is 
introduced. 


From  the  literature  on  the  principles  of  communication, 
Mathews  and  Aghvami  draw  general  bounds  on  the  BER.  These 
are  valid  for  any  type  of  signal  corruption  and  can  be 


written 


h(P1+P2) <PM(e) <P1+P2 


(122) 


P-^(e)  =  Pr  [  —  ao <^q  <_  0 ] 


(123) 


P2 (e)  =  Pr  [  -c 


><E  <j« ;  tan  — <_Q  £  °°J 


(124) 


For  linear  CPSK  systems  with  noise  as  the  only  form  of  signal 
corruption,  the  in  the  lower  bound  is  replaced  by 
The  authors  propose  that  this  is  also  true  for  nonlinear 
systems  thus  making  the  upper  bound,  P-^  +  P2,  a  good  estimate 
of  the  true  error  probability.  This  estimate  is  used  to 
develop  a  generalized  conditional  error  probability. 


PM(e|R,S) 


fe)  -er£  ( 


(125) 


where  IQ  and  Qq  are  I  and  Q  without  the  vc  and  vg  terms 
respectively.  By  averaging  over  all  possible  values  of  R 
and  6,  the  system  error  probability  is  obtained.  This  last 
step  is  relatively  uncomplicated  and  straightforward  as  the 
joint  pdf  of  the  amplitude  and  phase  of  a  known  sinusoid 
corrupted  by  a  narrowband  Gaussian  noise  process  [q(R,3)]  is 
well  known. 

Various  and  sundry  mathematical  manipulations  lead  to 
the  final  form  of  the  BER  expression  that  the  authors  choose 
to  use — albeit  a  rather  arbitrarily  chosen  one.  Picking 


expressions  for  f(>)  and  g(*)  that  are  used  m  other  works, 
Mathews  and  Agvami  proceed  to  compute  PM(e)'s  for  M  =  2,  4, 

8,  and  16  and  present  the  graphs  for  the  latter  two.  For 
the  binary  and  quaternary  cases,  it  is  found  that  the  results 
computed  are  in  exact  agreement  with  those  from  [18].  When 
the  signal  is  octernary  or  duo-octernary  and  the  proposed 
bound  is  applied  using  both  uplink  and  downlink  E^/N^'s, 
the  results  are  seen  to  be  "good"  approximations — less  than 


10%  deviation  from  the  true  PM(e)'s, 


The  authors  conclude 


that  their  method  is  acceptably  accurate,  much  simpler  than 
the  calculation  of  the  exact  BER,  and  therefore,  much  better 
than  the  analytical  approaches  whose  culmination  is  in  the 
work  of  Forsey  and  associates. 

The  most  recent  analytical  approach  to  a  wideband  non¬ 
linear  system  with  significant  uplink  noise  is  described  by 
Yao  and  Milstein  in  a  1982  paper  [20].  Contending  that  the 
final  test  statistic  cannot  be  adequately  approximated  as 
the  sum  of  2BT  independent  samples,  nor  by  a  majority  logic 
detector  utilizing  the  sgn ( • )  function,  nor  by  a  single 
sample,  as  all  of  the  previously  examined  work  has  done,  the 
authors  use  the  moments  of  the  system's  output  to  predict 
system  BER  performance.  With  the  same  general  approach  used 
as  that  found  in  the  work  of  Milstein  and  Davidson  [14],  a 
numerical  procedure  based  upon  the  use  of  moment  space  bounds 
and  introduced  in  earlier  work  by  these  authors  is  presented. 

Following  common  lines  of  system  and  signal  development, 
the  average  error  probability  of  the  system  is  seen  to  be 


PM(e)  *  E 


[q  ‘(ili  “i/^) 


(126) 


2m2 
a2  =  N  av 


(127) 
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where  is  a  sample  of  the  demodulated  input  to  the  final 
integrate  and  dump  detector  due  to  the  signal  component  of 
the  received  waveform.  Utilizing  the  moment  space  bounding 
technique  and  more  specifically,  the  isomorphism  theorem, 
two  moments  bounding  the  error  probability  are  selected. 

Through  a  process  not  well  understood  by  this  researcher, 
a  group  of  generalized  equations  are  developed  that  will, 
with  the  numerical  determination  of  two  sequences  of  constants, 
be  applicable  to  systems  with  arbitrary  nonlinearities. 

Yao  and  Milstein  proceed  to  present  numerical  results 
for  the  case  of  an  ideal  hard-limiter  as  the  system  non¬ 
linearity  for  various  values  of  uplink  and  downlink  Eb/NQ. 

3.1.3  Nonlinear  Systems  with  Uplink 

Noise  and  IS I 

The  last  of  the  analytical  methods  to  be  discussed  are 
those  that  are  used  on  more-or-less  realistic  systems.  That 
is  to  say,  those  approaches  that  have  been  applied  to  the 
basic  PSK  system  as  described  at  the  beginning  of  this 
chapter.  Unfortunately,  the  quantity  of  work  within  this 
classification  is  quite  small  comprising  only  four  papers — 
three  of  which  were  jointly  authored  by  N.  Ekanayake  and 
D.  P.  Taylor. 

In  their  earliest  work  [27]  (the  first  to  be  written, 
although  not  the  first  to  be  published) ,  Ekanayake  and  Taylor 
present  what  is,  essentially,  just  an  extension  of  the  work 
of  Jain  and  Blachman  [16]  to  M-ary  systems  with  ISI.  Like 
Jain  and  Blachman,  the  authors  assume  that  the  transponder 
nonlinearity  is  an  ideal  hard-limiter.  Additionally,  all  of 
the  ISI  is  assumed  to  arise  at  the  transmit  filter  while  the 
receive  filter  is  wideband  and  introduces  no  signal  corruption. 
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Following  the  standard  signal  development  up  to  the 
transponder  input,  the  authors  define  an  "effective"  pulse 
shape  using  the  quadrature  signal  components  such  that 

s  (t-kTs)  cos  v(t-kTs)  =  pk  cos  4>k  +  i^k  p^cos  4^  +  nc, 

(128) 

s(t-kTs)  sin  v(t-kTs)  =  pk  sin  $k  +  ^  p±sin  4^  +  ng . 

(129) 

This  gives  the  satellite  input  the  form 

Xp  =  k  sk  cos(u)0t  +  (130) 

which  leads  to  the  transponder  output  being  a  cosinusoid  of 
frequency  u  ,  phase  v,  ,  and  an  amplitude  of  unity.  In  the 
receiver,  the  signal  is  coherently  demodulated  and  sampled 
as  in  [18]  yielding  inphase  and  quadrature  signal  components 
that  are  written  as 

I  =  cos  vk  +  vc  (131) 

Q  =  sin  vk  +  vs 

for  any  k-th  time  period.  It  is  quite  easy  to  see  that  the 
determination  of  the  probabilistic  character  of  v>k  has  become 
the  focus  of  the  authors'  work. 


The  derivation  of  the  error  probability  as  it  proceeds 
is  fairly  straightforward.  From  equations  (131)  and  (132), 
the  expression  for  the  probability  of  a  correct  decision 
[QM(e)]  is  written  where,  of  course,  PM(e)  =  l-QM(e). 

QM(e)  is  taken  as  an  expectation  over  v  and,  after  some 
algebraic  manipulations,  is  found  to  be  very  similar  in  form 
to  that  given  in  [16].  In  fact,  for  BPSK  with  no  ISI,  the 
expression  Q2  (e)  is  seen  to  be  identical  to  that  derived  by 
Jain  and  Blachman.  At  this  point,  the  only  unknown  is  the 
pdf  of  v. 

Since  the  next  logical  step  is  the  calculation  of  q(’j), 
that  is  what  Ekanayake  and  Taylor  proceeded  to  do.  Using 
techniques  and  solutions  previously  introduced  by  Shimbo, 
et  al.,  in  [12],  an  expression  for  q(v)  is  given  and  the 
moments  required  for  the  calculation  of  the  BER  are  presented. 
Following  these  steps,  the  manipulations  required  to  achieve 
a  relatively  simplified  expression  for  PM(e)  are  performed. 

The  paper  is  concluded  with  several  numerical  examples 
utilizing  the  transmit  filter  used  in  [12]  and  a  single 
sample  per  symbol  receiver.  Various  graphs  are  presented  for 
BPSK  and  QPSK  systems  with  different  BTg  products.  Additionally, 
each  graph  shows  the  BER  for  the  case  when  no  ISI  is  present. 

For  large  values  of  the  effective  uplink  E^/N^  (>  lOdB)  it 
is  noted  that  the  calculation  of  PM(e)  becomes  quite  tedious 
and  an  asymptotic  approximation  is  given  so  that  computational 
time  may  be  reduced. 

A  final  note  must  be  made  about  the  work  just  descirbed: 

For  M  >  2,  the  authors  made  a  mistake  in  the  derivation  of 
q(v).  This  matter  is  acknowledged  and  corrected  in  an  errata 
announcement  published  in  the  "IEEE  Transactions  on  Information 
Theory,"  vol.  IT-27,  no.  1,  1/81,  pp  137-138.  Although  the 
correction  introduces  Hermitian  functions  not  found  in  the 
original  work,  the  basic  premises  are  found  to  be  quite  sound. 


In  their  work  immediately  following  that  presented  in 

[27],  Ekanayake  and  Taylor  present  a  method  of  predicting 

system  performance  that  is  a  step  backward  as  far  as  accuracy 

is  concerned.  Recognizing  the  computational  complexity  of 

the  method  presented  previously — one  that  requires  both  an 

exact  knowledge  of  the  transmit  filter  impulse  response  and 

the  evaluation  of  an  infinite  series — the  authors  present 

simple  upper  and  lower  bounds  on  the  BER  of  a  BPSK  system 

with  ISI  and  a  hard-limiting  transponder.  Since  the  hard- 

limiter  is  assumed  to  be  ideal  and  with  the  presumption  that 

no  ISI  arises  in  the  receive  filter,  these  bounds  are  very 

easy  to  apply  requiring  knowledge  only  of  the  peak  value  and 

variance  (c^  )  of  the  intersymbol  interference.  Estimates  of 

these  two  quantities  are  generally  easy  to  find  and,  in  fact, 

can  be  taken  from  some  of  the  curves  developed  in  [16], 

The  development  of  the  signal  is  pursued  in  much  the 

same  manner  as  has  been  seen  previously.  One  of  the  differences 

is  that  the  ISI  is  assumed  to  be  finite  and  to  exist  between 

the  limits  of  +1  and  -I  where  I  is  the  maximum  possible 

m  m  m 

amplitude  of  ISI  caused  by  the  N  pulses  preceding  and  following 
the  given  pulse.  Another  is  that  the  uplink  noise  and  the 
downlink  noise  are  both  given  the  general  form 

n(t)  =  n  ( t )  cos(u  t+$.  )-n  (t)  sin(cu  t+<j>,  )  (133) 

G  O  K.  S  O  iv 


thus  causing  the  combined  phase  at  the  transponder  input  to 
be  written  as 


\k(t)  =  arctan 


r 

ng 

V  *  P? 

i/k 


(134) 


Since  the  limiter  is  assumed  to  be  bandpass  and  ideal,  the 
output  is  just  a  cosinusoid  with  phase  (4>^  +  •  The  signal 


is  subsequently  corrupted  by  downlink  noise  and  then  coherently 
demodulated  and  sampled  at  the  receiver. 

With  y  defined  as  the  numerator  of  (134)  and  x  as  the 
denominator,  the  average  probability  of  error  is  written 


?2  (e)  =  Ex^[Js  erfc  (cos  \^//2  a^)  ] 


(135) 


and  evaluated  by  integrating  over  all  possible  values  of  the 
ISI  with  respect  to  its  distribution  function  at  the  sampling 
instant.  A  little  thought  leads  one  to  the  fact  that  the 
evaluation  of  (135)  is  accomplished  quite  easily  in  the  absence 
of  ISI  and  is  therefore  wholly  dependent  on  the  unknown 
distribution  of  the  ISI.  If  the  ISI  is  identified  by  x, 
and  the  error  probability  without  ISI  effects  is  denoted  by 
Pg(pk+a) ,  then  some  reflection  on  the  nature  of  (135)  leads 
to  the  bounds 


M 

h[P  (Pt+I  )  +  p  (Pt,-I  )]  -  p  (PtJ|  +  P  (P,  ) 

g  He  m  g  He  m  g  'He  1  I  gvtK 


1  p2(e)  ^Pg(Pk+ml)  +  Pg  (pk_ml)  1 


(136) 


where  is  the  first  absolute  moment  of  the  ISI.  Since  m^ 
itself  is  very  hard  to  specify,  the  bounds  introduced  for  it 
in  [8]  are  used.  Using  curves  developed  in  [16],  (136)  can 

be  computed  and  is  for  a  4th  order  Che’ yshev  filter  by  way 
of  example. 

The  authors  conclude  that  the  bounds  presented  are 
acceptably  tight  so  long  as  (E^/N^) u  is  greater  than 
( Eb/NQ )  d *  In  9enera1'  it  is  noted  that  the  lower  bound  is 
tighter  than  the  upper.  Further,  the  divergence  that  does 
occur  when  downlink  E^/N0  exceeds  that  of  the  uplink  is  seen 
to  approach  a  constant  value  with  increasing  (E^/N^^. 
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In  the  two  papers  just  described,  Ekanayake  and  Taylor 
present  two  ways  to  predict  a  satellite  communications 
system's  BER  performance.  The  first  allows  an  exact  solution 
whereas  the  second  is  a  presentation  of  acceptable  bounds 
designed  to  sidestep  the  computational  complexity  of  the 
exact  solution.  In  [29],  the  authors  present  the  third  paper 
on  their  work  in  this  field.  There  are- no  new  approaches, 
however,  as  the  reader  is  merely  given  a  rehash  of  work  that 
has  already  been  discussed  in  [28].  In  fact,  [29]  is  just 
an  explication  and  justification  of  the  work  presented  in 
the  authors'  previous  paper  and,  as  such,  need  not  be  further 
explained. 

3.2  Monte  Carlo-Based  Simulations 

True  Monte  Carlo  (MC)  simulation,  although  the  most 
accurate  predictor  of  a  communications  system's  performance, 
calls  for  an  unacceptably  large  amount  of  computational  time 
and  effort  to  produce  reliable  estimates  of  the  BER.  There¬ 
fore,  the  technique  is  not  used  in  its  pure  form.  It  is, 
nonetheless,  important  to  understand  the  ideas  and  procedures 
underlying  MC  simulation  as  a  great  deal  of  the  work  on 
satellite  communications  system  performance  prediction  is 
based  upon  these  very  same  concepts. 

The  estimation  of  the  probability  of  error  for  an  M-ary 
version  of  the  system  shown  in  Figure  1  utilizing  pure  MC 
techniques  proceeds  as  follows:  First,  one  generates 
sequences  of  sampled  values  of  D(t),  n(t),  and  v(t)  using 
random  number  generators  thus  allowing  the  use  of  deterministic 
models  of  the  various  system  components  and  elements.  The 
sequence  { }  passes  through  the  transmit  filter  becoming 
{ S ^ }  such  that 


where  one  can  easily  see  that  the  system  is  assumed  to  have 
a  memory  length  of  J  representing  the  ISI  on  the  uplink. 

J  is  selected  so  that  the  vast  majority  of  each  pulses's 
energy  is  accounted  for.  Following  the  by  now  familiar 
progression  of  the  signal  through  the  system,  one  arrives 
at  the  input  to  the  receive  filter  {U^K  The  k-th  output 
of  this  filter  is  expressed  as 

Yk  =  hr(Uk'  Uk-1 '  *  *  *  '  S-R-l*  (138) 

where  R  fulfills  the  same  function  for  the  downlink  as  J 
does  for  the  uplink.  R  and  J  may  be,  but  are  not  necessarily, 
the  same. 

At  this  point,  {Y^}  is  sampled  at  the  system  bit  rate 
and  compared  to  the  decision  threshold  or  thresholds  to 
determine  the  output  bit.  The  system  error  counts  are  obtained 
by  comparing  the  output  of  the  receiver,  {b£},  with  the  known 
input  bit  sequence  (b^).  In  order  to  estimate  the  error 
probability  by  this  direct  counting  method,  the  number  of 
samples  required,  N,  is  given  by 

N  >  -~-i -  (139) 

s  ve) 

where  e  is  the  normalized  error  of  the  estimated  P^(e).  That 
is, 

(140) 


standard  deviation  of  P„(e) 


It  is  quite  easy  to  see  that  the  minimum  number  of  bits 
required  for  accuracy,  N  =  10/PM(e),  yields  a  normalized 
error  of  e  =  . 32  %  j  —  an  error  that  may  be  too  high  in  many 
cases  thus  requiring  even  larger  sample  sizes.  The  uselessness 
of  pure  Monte  Carlo  simulation  is  obvious. 

In  light  of  this  problem,  the  objective  of  all  of  the 
simulation-based  performance  prediction  tools  is  to  drastically 
reduce  the  sample  size  requirement  without  losing  any  of  the 
flexibility  inherent  in  Monte  Carlo  simulations.  The  pursuit 
of  this  goal  has  resulted  in  the  appearance  of,  roughly,  three 
categories  into  which  the  vast  majority  of  simulation  tools 
may  be  classified.  For  the  following  descriptions,  these 
categories  shall  be  known  as  1)  the  Combinatory/Short  Sequence 
or  Semi-Analytical  approach,  2)  the  Tail  Probability/Extreme- 
Value  Theory  approaches,  and  3)  the  Modified  Monte  Carlo  or 
Importance  Sampling  approach. 


3.2.1  The  Semi-Analytical  Approach 

In  the  design  of  a  simulation-based  performance  prediction 
tool,  the  foremost  goal  is  the  reduction  of  the  required 
sample  size.  What  makes  this  goal  difficult  to  realize  is 
the  concurrent  requirement  that  the  accuracy  and  flexibility 
of  a  true  Monte  Carlo  simulation  be  retained.  The  primary 
cause  of  the  large  sample  size  requirement  is  the  necessity 
of  modeling  the  uplink  and  downlink  noise  using  sampled  random 
processes  the  elimination  of  which,  therefore,  would  seem  to 
be  the  answer  to  the  problem.  Any  analytical  representation 
of  noise  will,  of  course,  cut  down  on  the  flexibility  of  the 
tool;  however,  most  noise  that  is  actually  encountered  can 
be  represented  by  one  of  several  well  known  analytical  noise 
models  thereby  bypassing  this  objection.  Since  the  downlink 


noise  encounters  only  linear  system  elements  there  is  no 
problem  in  using  an  analytical  model  for  its  representation. 

The  uplink  noise  must  pass  through  a  nonlinear  transponder 
however,  and  one  thereby  encounters  trade-offs  in  simulation 
time  and  accuracy  when  an  analytical  treatment  is  to  be 
utilized.  Replacing  the  random  process  representations  of 
the  uplink  and  downlink  noise  with  some  form  of  analytical 
modeling  and/or  treatment  while  attempting  to  retain  the 
flexibility  and  accuracy  of  MC  simulation  yields  what  may  be 
designated  as  the  semi-analytical  method. 

Although  several  slightly  different  approaches  were 
reviewed,  all  were  essentially  the  same  and  the  representative 
papers  on  the  subject  [31] -[36]  may  all  be  treated  and  described 
as  one  approach.  That  approach  is  to  completely  remove  all 
random  processes  from  within  the  framework  of  the  simulation. 

The  input  sequence  is  handled  by  using  several  short  sequences 
where  the  system  BER  becomes  the  average  of  the  BERs  calculated 
for  each  sequence  so  that  the  ISI  effects  are  still  accounted 
for.  In  each  of  the  performance  prediction  tools  described 
in  [31] -[36]  all  of  the  noise  in  the  system  is  presumed  to 
arise  on  the  downlink.  Some  use  the  rationale  that  since  the 
uplink  S/N  is  typically  much  larger  than  that  of  the  downlink 
then  the  effect  the  uplink  noise  has  on  the  probability  of 
error  is  so  slight  as  to  be  insignificant  and  can  be  safely 
ignored.  Others  assume  that  the  uplink  noise  may  be  significant 
and  compute  an  "equivalent"  downlink  noise  power  using  some 
simple  algorithm  in  a  manner  vaguely  similar  to  the  work  of 
Jain  and  Blachman  [16]  or  Lyons  [15].  In  short,  the  semi- 
analytical  method  is  comprised  of  the  retention  of  the  "building 
block"  analytical  models  of  Monte  Carlo  simulation,  the  use 
of  several  short  sequences  as  input  to  replace  that  random  or 
partially  random  signal,  and  the  analytical  treatment  of  the 
thermal  noise  as  if  it  were  all  downlink  and  thus  easily  handled 
in  a  deterministic  manner.- 


3.2.2 


Probability  Approaches 

The  preceding  section  described  the  commonly  used 
approach  to  the  prediction  of  a  communications  system*  s 
performance  that  entails  the  alteration  of  a  system’s 
simulation  model  ,so  that  the  nonlinear  transformation  of  the 
uplink  noise  can  be  handled  without  resorting  to  a  true 
Monte  Carlo  simulation.  An  alternative  or  some  alternatives 
to  that  approach  is  to  assume  that  a  limited  amount  of  Monte 
Carlo  data  is  to  be  taken  and  to  use  these  data  to  predict 
the  system's  BER.  Two  such  alternatives  encompass  the 
application  of  the  theory  of  the  asymptotic  distributions 
of  extreme  values  [37]  —  [42] . 

Of  the  two  alternatives  mentioned  above,  the  earliest 
to  be  proposed  was  the  direct  and  generalized  application 
of  the  statistical  theory  of  extreme  values  to  communications 
system  performance  prediction  [37]— [39] ,  [42].  Suggested 

for  and  described  in  the  context  of  a  binary  system  where 
one  is  making  just  one  threshold-crossing  decision,  it  is 
applicable  to  more  complex  systems  with  a  corresponding 
increase  in  complexity.  With  the  primary  assumption  that 
the  output  and/or  degradation  of  the  system  output  can  be 
represented  as  a  random  process  possessing  a  distribution 
function  of  exponential  type  on  the  right,  the  application  of 
extreme-value  theory  is  described  in  the  next  few  paragraphs. 

Let  one  assume  for  explicative  reasons  that  the  transmitted 
symbol  (or  symbols)  is  (are)  0.  Further,  let  {X^)n  represent 
a  sequence  of  n  independent  samples  drawn  from  the  distribution 
Fx(x)  of  the  output.  If  the  random  variable  Z  is  defined  as 


Z  =  maximum  (X.,  ,  .  .  .  , 


XJ 


(141) 


»l»  I 


then  the  distribution  function  Fz(z)  is  said  to  take  the 

asymptotic  form 

« 

lim  Fz (z)  =  lim  Fxn(z)  =  exp (-exp [-a ( z-u) ] ) .  (142) 

n  -*■  09  n  -*■  ® 

this  lends  itself  to  a  linearization  yielding 

-ln(-ln[Fz (z) ] )  =  a(z-u).  (143) 


a  and  u  in  the  previous  two  equations  are  unknown  parameters 
of  the  asymptotic  distribution  that  must  be  estimated  from 
the  Monte  Carlo  data  that  is  presumed  to  exist.  Since  the 
quantity  -ln(-ln[Fz (z) ] ) ,  and  therefore,  the  error  probability 
for  a  given  value  of  z,  is  now  linear  in  z  we  can  use  this 
theory  to  estimate  a  system's  BER  performance. 

Extreme-value  theory  performance  prediction  begins  by 
dividing  the  Monte  Carlo  data  into  N  groups  of  m  samples. 
Although  there  are  no  optimum  procedures  specifying  the 
selection  of  N  and  m,  the  theory's  requirements  call  for  m 
to  be  large  enough  that  the  extreme  samples  have  the  asymptotic 
distribution  of  (143)  and  for  N  to  be  large  enough  that  the 
extrema  from  each  group  of  m  samples  can  be  used  to  acquire 
unbiased  estimates  of  the  parameters  a  and  u.  Following  the 
division  of  the  data,  the  estimates  of  a  and  u  (a'  and  u") 
are  determined  as  follows: 


1.  Select  the  largest  sample  from  each  of  the  N 
groups  X(]_),  X(2)'*-.X(N)  and  order  them  so  that 
x(l)  1  x(2)  1  x(3) ... 

2.  Estimate  the  value  of  Fz(x^)  as 


N+l-k 


N 


FX(x,J 


\T 


k  =  1 ,  2 , . . . 


(144) 


3.  Letting  y^  =  -ln(-ln[F2 (xk) ] ) ,  Plot  a  "best"  linear 
fit  to  the  data  (x^y^)  and  draw  the  estimates  for 
a  and  u  from  the  resulting  coefficients. 

At  this  point  one  may  obtain  the  estimated  value  of  the 
distribution  at  the  decision  threshold  level  T  as 


F'(T)  =  exp(-expt-a>  (T-u')  ]  ) . 


(145) 


A  diagram  of  this  is  shown  in  Figure  8.  Now  we  know  that 


VT)  =  FX‘“(T) 


«  Cl-P2(e)]1 


(146) 


%  1  -  mP2(e) 


and  the  estimate  for  the  probability  of  error  is  therefore 


P2(e)  = 


1  -  F'(T) 


(147) 


an  estimate  which,  furthermore,  is  consistent  and  asymptoti¬ 
cally  unbiased. 

Study  and  analysis  of  the  accuracy  and  of  the  effectiveness 
of  this  method  are  given,  usually  assuming  a  maximum- 
likelihood  receiver.  The  accuracy  is  deemed  quite  acceptable 
and  the  effectiveness,  that  is,  the  sample  size  reduction 
factor  from  that  required  for  a  true  Monte  Carlo  simulation, 
is  shown  to  be  8  for  an  error  rate  of  10  ^  with  a  confidence 
level  of  90%  and  a  minimum  resolution  of  less  than  one  order 
of  magnitude.  This  factor  is,  of  course,  expected  to  be  much 
larger  for  higher  BER's  and/or  confidence  levels. 


Figure  8  -  Extreme-Value  Theory  Sample  Diagram 
(for  N=4) 


Following  the  basic  presentation  of  this  given  in  [37], 
[38]  and  [39]  present  generalizations  of  the  theory  of  extreme 
value  distributions  in  what  are  essentially  efforts  to  prove 
the  validity  of  the  application  presented  by  Posner  while 
pointing  out  some  of  the  mistakes  and  misinterpretations 
commonly  made.  It  is  noted  that  the  main  problem  encountered 
is  that  the  convergent  sequence  in  (142)  is  frequently  assumed 
to  apply  when,  in  fact,  there  is  no  convergence.  This,  of 
course,  leads  to  inaccurate  and  bad  estimation  of  system 
error  probabilities.  Letting  w  =  a(z-u)  and  assuming  the 
existence  of  some  sequence  {d.  },  Weinstein  [38]  rewrites 
equation  (142)  as 


lim 
n  -*■  o 


Fz(z)  =  lim 


V11 


n  -*■  ® 


(•’  *  *) 


(148) 


=  exp [-exp (-w) ] 


where  V(*)  is  probability  distribution  and  v  >  0.  With  the 

introduction  of  this  form,  he  proceeds  to  define  the  necessary 

and  sufficient  conditions  for  the  convergence  of  (148).  v  has 

been  introduced  primarily  as  a  generalizing  factor  in  the 

equation  that  more  adequately  shows  the  true  nature  of  the 

convergence  to  the  extreme-value  asymptote.  He  proceeds  to 
* 

show  how  (148)  leads  to  a  more  consistent  and  accurate  BER 
estimator.  In  [39],  Jeruchim  elucidates  the  applications 
of  this  generalized  estimator  and  its  potential  pitfalls. 

The  second  alternative  to  be  based  upon  the  theory  of 
the  asymptotic  distribution  of  extreme  values  is  presented 
by  Weinstein  and  follows,  albeit  sc  .iewhat  indirectly,  from 
the  generalization  he  presented  in  [38].  Noting  that  in  some 
cases  substantial  inaccuracies  are  inherent  in  the  use  of  an 


asymptotic  distribution  instead  of  the  actual  extreme-value 
distribution,  he  derives  a  BER  estimator  that  is  very 
specifically  directed  towards  communications  systems. 

Weinstein  begins  with  the  assumption  that  there  is  a 
general  class  of  probability  distributions  that  most  communi¬ 
cations  system  degradation  may  be  described  by.  With  the 
requirement  that  this  class  of  distributions  be  linearized 
by  a  known  transformation  into  a  simple  function  of  the 
argument,  he  describes  these  distributions  by  the  probability 
density  function 


qa,a  <X)  =  — (W  exP(-<x/°>a1 


(149) 


where  r(-)  is  the  Gamma  function.  The  random  variables 
described  by  (149)  are  called  the  higher  exponential  class 
of  random  variables. 

Considering,  once  again,  a  binary  system,  the  probability 
of  error  can  be  written  as 


P2(e)  =  f  %,a(x>  dx 


(150) 


where  T  is  the  threshold  and  a  0  is  presumed  to  have  been 
transmitted.  The  well  known  asymptotic  expansion  of  this 
equation  leads  to 


P2(e)  £  exp  (-Ta) 


T  >>  1. 


(151) 


which  is,  of  course,  linearized  by  a  double  logarithmic 
transformation  to  the  form 


In (-In [P- (e [ 0)  ])  £  a  ln(T) • 


(152) 


Taking  this  linearization,  Weinstein  uses  an  extrapolation 
technique  similar  to  that  described  earlier  in  this  section 
to  acquire  an  estimate  of  the  error  probability.  The 
procedure  is  as  follows: 

1.  Select  three  extrapolation  thresholds,  x^,  X2 ,  and 
X3  such  that  a  small  portion  of  N  Monte  Carlo  data 


observations  exceeds  x- 


The  three  thresholds 


should  be  selected  to  be  un.-1'  formly  spaced  in 
logarithmic  scale. 

2.  Calculate  the  estimated  probability  of  error, 
p2  (e  |  x^) ,  for  each  of  the  thresholds  as 


p2(e|xk)  -  N 


(153) 


where  n^  is  the  number  of  observations  that 
exceed  x^. 

Plot  a  "best"  linear  fit  to  the  data  [x^,  PfCejx^)] 
and  draw  the  estimate  of  a  from  the  coefficients. 

At  this  point  the  estimated  probability  of  error 
for  the  threshold  T  can  be  found  from 


Pi  (e  I  T)  =  exp  (-Tc 


(154) 


Examinations  of  this  estimator  have  shown  that  with 
proper  selection  of  x^,  x 2,  and  x^,  it  is  acceptably  accurate, 
insistent  and  relatively  unbiased.  There  are  no  procedures 
for  the  selection  of  the  optimum  values  of  the  x^,  and  they 
are  usually  chosen  as  a  compromise  between  an  excessive 
sampling  variance  when  they  are  too  large  and  an  unacceptable 
amount  of  estimation  bias  when  they  are  too  small. 

The  sample  size  reduction  factor  for  this  "tail  probability" 
approach  is  shown  to  be  roughly  3.8  when  a  =  1,  18  when 
a  =  2,  and  54  when  a  =  3.  In  the  binary  system  case,  an 
additional  factor  of  2  applies  for  all  values  of  a.  These 
numbers  are  derived  for  a  probability  of  error  of  10 
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The  Importance  Sampling  Approach 


The  last  of  the  three  simulation  approaches  is  what  is 
known  as  the  importance  sampling  or  modified  Monte  Carlo 
approach  [43]— [46] -  Like  the  extreme-value  theory  approaches 
importance  sampling,  as  applied  to  performance  prediction, 
is  based  upon  true  Monte  Carlo  simulation  methods.  That  is, 
the  communications  system  BER  is  estimated  from  a  set  of 
Monte  Carlo  observations  rather  than  using  the  analytical 
techniques  described  in  Section  3.2.1. 

Importance  sampling  is  applied  to  the  same  basic  system 
that  has 'been  discussed  throughout  this  chapter  and  in  the 
preceding  chapters.  The  formulation  of  the  problem,  however, 
can  be  stated  in  a  somewhat  different  manner.  To  begin, 
one  assumes  that  the  input-output  relationship  of  the  communi¬ 
cations  system  in  question  is  known.  In  addition,  the 
statistics  of  the  inputs  to  the  system  are  assumed  to  be 
known. 

Our  desire  is  to  use  this  knowledge  to  estimate  the 
probability  density  function  of  the  output  from  which  we  can 
derive  an  estimate  of  PM (e ) .  If  the  system  is  once  again 
assumed  to  be  binary,  then  the  error  estimate  can  again  be 
written  as 


p2(e) 


(155) 


where  <  is  the  decision  threshold  and  t(.  is  the  estimated 
pdf  of  the  output  sequence  {Y^}.  It  has  already  been  presumed 
that  the  error  estimate  will  come  from  a  group  of  N  Monte 
Carlo  observations.  This  means  that 


where  is  the  number  of  the  output  observations  that  exceed 
<.  The  number  of  output  samples  required  to  estimate  the 
error  probability  via  Monte  Carlo  simulation  is,  however, 
known  to  be  given  by 


N  >  l/e2PM(e) 


(157) 


with  e  being  the  normalized  error  of  the  estimate  P2(e) 

(for  £  =  1//10  this  yields  N  >  10/PM(e)  mentioned  earlier) . 
Since  (157)  requires  that  N  be  unacceptably  large  for  an 
acceptably  small  e  given  typical  communications  systems  error 
probabilities ,  the  problem  becomes  one  of  estimating  f£  to 
minimize  e  for  a  reasonable  and  acceptable  number  N  of 
observations  of  the  system  output. 

In  the  application  of  importance  sampling  to  the  problem 
stated  above,  the  principle  is  to  reduce  the  sample  size 
requirement  by  modifying  the  statistical  properties  of  the 
systems  input  sequence.  This  causes  a  significantly  higher 
than  normal  number  of  errors  to  occur  in  the  course  of  the 
Monte  Carlo  simulation.  This  modification  must  be  compensated 
of  course,  by  weighting  the  error  estimation  of  the  true  BER. 

The  quantity  that  equation  (155)  is  estimating  is  written 


P9(e)  -  /  fv  (y)  dy. 

0  x 

It  can  be  rewritten  as 

as 

P2(e)  -  /  d(y)  fY  (y)  dy 


(158) 


(159) 


where  we  have  defined  d(y)  =  1  if  y  exceeds  the  threshold 
6  and  0  otherwise.  In  importance  sampling,  a  modified  pdf 


f*Y(y)  is  introduced  that  is  superior  for  sampling  purposes. 
Equation  (159)  can  therefore  be  rewritten 


P2(e)  mJ,mdtY)  r^Ty T  fVy)  dy 

* y  d*(y)f*Y(y)  dy. 


(160) 


The  modified  output  pdf  is  derived  by  biasing  the  pdf  of  the 
input  and  (156)  now  becomes 


P^(e)  =  i 


N 

E  d* 
i=l 


(yi> 


(161) 


with  y^  the  i-th  element  of  the  output  sequence  {Y^}  that 
is  produced  by  using  the  biased  input  sequence. 

In  the  implementation  of  the  modified  Monte  Carlo 
technique  one  follows  essentially  the  same  procedures  as  in 
true  Monte  Carlo  simulation.  There  are,  however,  some  rather 
obvious  differences.  From  (160)  it  is  seen  that  the  output 
samples  will  have  a  bias  that  is  given  by  the  ratio  f*y(y)/ 
fY(y).  It  is  therefore  obvious  that  the  input  samples  will 
have  a  probability  bias  B(x)  which  is  the  ratio  f*x(x)/fx(x) . 

The  weight  that  must  be  given  to  the  effect  of  the  sample, 

W(x)  is,  therefore,  the  inverse  of  the  bias.  Since  (Y^)  is 
determined  from  {X^},  the  probability  that  the  i-th  sample 
of  Y  *  y.  is  also  assumed  to  have  the  bias  B(x^)  and  thus 
the  weight  W(xi).  The  probability  estimate  determined  using 
equation  (156)  must  therefore  be  weighted  by  the  average  of 
the  weights  of  all  of  the  n^  samples  that  are  observed  to  exceed 
the  threshold  [as  shown  by  equation  (161)].  In  other  words. 
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E  W (x. • ) 
nj  i=l  « 


(162) 


If  one  now  assumes  that  there  exist  N  observations  of 
the  output  of  a  system  having  an  input  whose  pdf  has  been 
biased  by  B(x),  then  equations  (156)  and  (161)  yield  an 
estimate  of  the  probability  of  error  as 


*j<«> 


wavj  VN 


■  w<xij> 


(163) 


which  can  be  shown  to  be  an  unbiased  estimator:  The  advantage 
of  using  this  modified  estimator  is  that  for  a  particular 
sample  size  the  variance  from  equation  (163)  will  be  smaller 
than  that  from  (156)  if  B(x)  is  chosen  to  be  >1  for  every 
value  of  x  e  I.  =  {xjy  *  g(x)e[y  >  8]}.  This  variance 
reduction  may  be  computed  as  follows. 

The  variance  of  the  estimator  of  equation  (156)  is  given 
by 

2  P2(e) [1-P' (e) ] 


*  P^eJ/N 


(164) 


f x (x)  dx* 


where  N„  is,  of  course,  the  number  of  observations  taken  for 
c 

the  application  of  the  Monte  Carlo  estimator  with  our  unbiased 
input.  The  variance  of  the  modified  estimator  of  (163)  is 
given  by 


ere*? 
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is2  4  B7-  ft.  w(x)  fx(x)  dx 


(165) 


where  W(x)  is  as  previously  defined  and  N^g  is  the  sample 

size  from  a  simulation  with  a  biased  input.  For  the  case  of 

N .  »  N  ,  the  variance  reduction  factor  becomes 

is  c 


fz  W(x)  fx(x)  dx 

_ 1 _ 

Si.  fX(x)  dx 


(166) 


Alternately,  if  it  is  assumed  that  the  normalized  error  as 
given  by  (157) ,  rather  than  the  sample  size,  is  held  constant, 
i.e.,  e  «  e.  ,  then  the  sample  size  reduction  factor  is 
given  by 


r  *  VNis 
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(167) 


From  the  previous  paragraph  it  may  easily  be  seen  that 
a  system's  BER  may  be  estimated  very  accurately  from  a  small 
group  of  observations  provided  that  the  input  pdf  has  been 
biased  such  that  those  portions  of  the  output  pdf  from  which 
errors  are  drawn  are  very  heavily  biased  [B(x)  >>  1]. 
Unfortunately,  causing  some  portions  of  the  output  pdf  to  be 
heavily  emphasized  causes  the  rest  to  be  heavily  de-emphasized. 
This  means  that  we  cannot  accurately  estimate  the  probabilities 
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of  error  that  are  drawn  from  the  de-emphasized  portions  of 
the  output  pdf;  however,  we  are  uninterested  in  those 
probabilities,  and  thus  this  is  not  a  problem. 

A  modified  Monte  Carlo  technique  utilizing  importance 
sampling  has  now  been  described.  The  implementation  can, 
however,  be  a  very  complicated  matter.  In  a  complex  system 
one  does  not  necessarily  know  beforehand  what  regions  of  the 
input  pdf  must  be  biased  to  produce  the  desired  bias  of  the 
output  pdf.  Fortunately,  in  most  communications  systems  the 
causes  of  performance  degradation  frequently  have  Gaussian 
probability  density  functions.  Therefore,  one  need  only 
bias  the  tails  of  the  input  pdf's  to  emphasize  the  desired 
portions  of  the  output  pdf's.  Methods  of  doing  this  and  the 
derivations  of  an  optimum  biasing  function  are  given  in 
[43]- [46] . 

In  addition  to  the  application  of  importance  sampling 
that  is  mentioned  above,  [43] -[46]  describe  the  extension 
of  the  technique  to  the  case  of  M  independent  system  inputs. 

This  is  shown  to  be  rather  simple  and  very  straightforward. 
Experimental  results  on  the  application  of  this  modified 
Monte  Carlo  technique  verify  its  usefulness  as  a  performance 
prediction  tool.  In  a  system  having  no  ISI,  the  sample  size 
reduction  factor  varies  from  45  -  25,500  depending  upon  the 
desired  probability  of  error  to  be  estimated.  With  ISI 
extending  across  five  symbols,  the  reduction  factor  is  shown 
to  be  on  the  order  of  7  -  1000.  Overall,  the  modified  Monte 
Carlo  approach  is  shown  to  meet  the  previously  stated 
requirements  with  its  advantage  over  true  Monte  Carlo  simulation 
increasing  as  the  desired  estimation  error  or  estimated 
P„(e)  are  decreased. 


4.0  PROPOSAL  AND  EXAMINATION  OF  A  PERFORMANCE 

PREDICTION  TOOL 

In  the  previous  chapter,  many  different  ways  to  predict 
the  performance  of  a  satellite  communications  system  were 
examined.  These  approaches  were  divided  into  two  main 
categories — the  analytical  methods  and  the  simulation-based 
methods — and  comprise  a  very  large  number  of  assumptions, 
mathematical  tricks  and  techniques,  and  unique  ways  of 
looking  at  the  " typical M  system.  Each  of  these  tools  were 
examined  to  see  if  they  satisfied  the  requirements  for  a 
performance  prediction  tool  as  stated  in  Chapters  1  &  2. 

Some  of  them  failed  because  they  did  not  address  the  effect 
of  a  nonlinearity  on  a  system's  performance.  Others  ignored 
the  bandlimiting  character  of  most  real  systems  or  else 
assumed  that  there  was  no  noise  added  prior  to  the  nonlinearity. 
Many  of  the  approaches  looked  at  in  Chapter  3  sacrificed 
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flexibility  for  accuracy  (and  vice  versa)  while  still  others 
were  so  computationally  complex  as  to  be  unfeasible.  Most 
failed  to  meet  the  requirements  for  combinations  of  these 
reasons. 

In  the  last  two  chapters  the  advantages  of  the  simulation 
based  approaches  over  the  analytical  approaches  have  been 
stated  time  and  again.  Although  this  may  not  be  true  for 
all  possible  simulation-based  or  analytical  methods,  it  is 
true  in  general  with  regards  to  the  research  that  has  been 
conducted  to  date.  Therefore,  all  further  discussions  shall 
center  entirely  upon  the  simulation-based  category  of  performance 
predictors. 

Of  the  three  subcategories  comprising  the  simulation- 
based  performance  prediction  tools,  each  has  its  advantages 
and  disadvantages,  yet  none  of  the  particular  methods  described 


in  the  previous  chapter  has  all  of  the  desired  properties 
of  the  ideal  prediction  tool.  The  semi-analytical  approach 
is  what  most  people  use  presently  and  thus  it  offers  little 
or  no  room  for  improvement.  Extreme-value  theory  is,  however, 
very  interesting  as  is  the  importance  sampling  approach. 

These  latter  two  techniques  offer  sample  size  reduction 
factors  ranging  from  3.8  up  to  many  thousands  while  retaining 
most  of  the  advantages  of  true  Monte  Carlo  simulation.  In 
fact,  the  only  significant  drawbacks  these  methods  have  are 
that  they  still  don't  reduce  the  sample  size  enough  for  a 
given  estimation  variance.  All  in  all,  these  methods  would 
seem  to  merit  further  examination. 


4 . 1  The  Proposal 

As  a  review,  let  us  state  briefly  the  basic  mechanisms 
by  which  Weinstein  and  then  Shanmugam  and  Balaban  achieve 
their  sample  size  reductions.  Weinstein's  method  is  one  of 
"graphic  extrapolation"  wherein  direct  counting  of  errors  is 
used  to  determine  three,  manageably  large  BER's  for  three, 
evenly  spaced  error  thresholds  x^  that  are  smaller  than  the 
threshold  of  interest  xR.  The  selection  of  the  x^'s  is 
governed  by  a  trade-off  between  the  sample  size  required  to 
accurately  estimate  the  BER  at  the  largest  x^,  x^,  and  the 
desire  to  have  x3  as  close  as  possible  to  xR  so  as  to  reduce 
the  variance  of  the  estimator.  This  trade-off  is  sample- 
size  limited  to  the  selection  of  an  x,  that  yields  a  PM(e) 
in  the  range  of  10~  to  10“  .  On  the  other  hand,  the  Modified 
Monte-Carlo  method  proposed  by  Shanmugam  and  Balaban  is  not 
strictly  sample-size  limited  at  all.  Their  method  requires 
that  one  have  a  prior  knowledge  of  the  statistical  character 
of  the  uplink  noise  especially  with  regards  to  that  portion 
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of  the  uplink  noise  pdf  that  directly  influences  the  number 

of  errors  to  be  seen  for  a  given  threshold.  A  little  thought 

will  indicate  that  a  reasonably  good  idea  of  the  statistical 

nature  of  the  output  signal  is  also  a  necessity.  With  the 

required  statistical  knowledge  in  hand,  one  biases  the  pdf 

of  the  uplink  noise  so  that  more  errors  occur  for  a  given 

threshold,  counts  the  errors  for  a  time  sufficiently  long  to 

yield  an  accurate  and  precise  estimate  and  then  removes  from 

this  estimate  the  bias  that  was  introduced  to  the  uplink  noise. 

The  sample  size  is  again  limited  to  that  needed  to  determine 
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a  BER  in  the  range  of  10  to  10  yet  the  determining  factor 
in  the  variance  of  the  estimator  is  more  the  accuracy  of  the 
biasing  applied  to  the  uplink  noise  statistics  than  in  how 
much  smaller  the  estimated  BER  is  from  the  counted  PM(e) 
used  in  the  estimation. 

With  the  mechanisms  of  these  two  methods  in  mind,  it  is 
time  to  directly  address  the  driving  force  behind  this  study: 
The  discovery  and/or  development  of  a  performance  prediction 
tool  for  satellite  communications  systems  that  more  closely 
approximates  the  "ideal"  performance  prediction  tool  described 
in  the  first  chapter  than  do  the  analytical  and  simulation- 
based  predictors  in  use  at  present.  Each  of  the  methods 
that  have  been  examined  has  been  shown  to  have  significant 
drawbacks.  With  the  exceptions  of  Weinstein's  method  and  the 
importance  sampling  approach  of  Shanmugam  and  Balaban,  each 
of  the  performance  prediction  tools  examined  has  been  shown 
to  be  unacceptable  for  one  reason  or  another.  To  be  thorough, 
in  fact,  not  even  the  two  exceptions  are  without  their 
problems.  In  short,  something  different  needs  to  be  examined 
and,  hopefully,  developed  to  fruition. 

A  different,  although  not  entirely  new  approach  to  the 
performance  prediction  tool  problem  will  be  offered  at  this 
time.  When  one  considers  the  two  approaches  that  were 
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previously  identified  as  being  the  best  of  those  that  are 
now  available  for  use,  it  is  easy  to  see  that  they  are  rather 
complimentary  in  nature.  That  is,  Shanmugam  and  Balaban's 
method  alters  the  statistics  of  the  uplink  noise  in  such  a 
way  that  either  the  number  of  samples  that  need  to  be  looked 
at  to  estimate  a  given  BER,  or  the  variance  of  an  estimate 
of  that  BER,  are  reduced  from  pure  Monte  Carlo  requirements 
whereas  Weinstein's  extrapolation  method  merely  lowers  the 
number  of  errors  that  need  to  be  counted  to  estimate  that 
given  BER.  It  should  be  a  fairly  straightforward  task  to 
take  the  two  methods  and  meld  them  into  a  new  tool  that  is 
based  upon  the  latter  and  uses  the  former  to  reduce  either 
the  sample  size  requirement  or  the  variance  of  the  estimator. 
In  short,  this  proposal  is  to  combine  the  two  best  approaches 
into  one  with  the  hope  of  counteracting  the  disadvantages  of 
each  while  enhancing  the  advantages. 

The  specifics  of  the  proposed  approach,  hereafter  to  be 
known  as  the  Modified  Extrapolation  Approach,  are  really 
quite  simple  and  are  as  follows: 


1.  One  begins  by  determining  the  expected  order  of 
magnitude  of  pM<e) . 

• 

2.  Define  a  maximum  number  of  samples,  Nme,  that  are 
to  be  taken  for  error  counting  purposes.  This 
will  be  of  the  available  computational  time  as  well 
as  the  minimum  number  of  errors  that  must  be  counted 
to  yield  a  reasonable  estimate  of  Pj^CelK)  where 

K  is  the  decision  threshold  (to  keep  the  variance 
small,  it  is  desirable  to  handle  as  close  to  the 
same  number  of  samples  as  required  by  Monte  Carlo 
estimation) . 

3.  Knowing  Nme,  the  largest  of  the  Weinstein  extrapolation 
thresholds,  *3me,  can  be  specified.  This  will 

depend  on  whether  the  modified  Monte  Carlo  method 
is  to  be  applied  to  further  reduce  the  sample  size 
or  to  lower  the  estimation  variance. 
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4.  From  X3me  determine  the  other  two  extrapolation 
thresholds. 

5.  Apply  the  bias  to  the  input  noise  pdf  and  start 
counting  errors. 


Evaluation  of  the  Proposal 


The  first  order  of  business  in  evaluating  the  proposal 
of  section  4.1  is  the  determination  of  criteria  by  which 
performance  is  to  be  judged.  Since  the  implementation  of 
the  Modified  Extrapolation  approach  can  be  broken  down  into 
the  separate  implementations  of  its  two  component  approaches, 
it  may  reasonably  be  assumed  that  the  selections  of  the 
extrapolation  thresholds,  the  bias,  B(x),  etc.,  can  be  based 
on  the  criteria  presented  by  Weinstein  or  Shanmugam  and 
Balaban  without  having  any  adverse  effect  on  the  performance. 

It  has  already  been  established  that  performance  prediction 
tools  are  judged  on  the  basis  of  accuracy,  flexibility,  and 
speed  of  execution.  By  the  nature  of  the  two  approaches  that 
are  being  combined  it  is  easily  seen  that  Modified  Extrapolation 
is  as  flexible  as  anything  that  now  exists  short  of  the  pure 
Monte  Carld)  technique.  We  are,  therefore,  left  to  consider 
the  accuracy  (variance)  of  the  estimate  and  the  speed  with 
which  that  estimate  is  made  (or,  stated  another  way,  the 
sample  size  reduction  that  is  achieved) . 

A  little  thought  will  suffice  to  convince  one  that  neither 
an  analysis  of  the  variance  (or  variance  reduction)  or  the 
sample  size  reduction  will  be  independent  of  the  other.  In 
order  to  study  the  effect  that  the  proposed  Modified  Extra¬ 
polation  method  has  on  the  variance  of  the  estimate  relative 
to  either  or  both  of  its  component  methods,  the  sample  sizes 
must  be  the  same,  and  vice  versa.  In  other  words,  testing 
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one  is  effectively  the  same  as  testing  the  other  given  that 
an  equal  importance  is  placed  on  the  two.  This  assumption 
has  been  and  will  continue  to  be  made.  Since  the  previous 
discussions  in  this  paper  have  tended  to  center  on  the  sample 
size  requirement  instead  of  the  variance,  the  following 
evaluation  of  the  Modified  Extrapolation  technique  will  do 
the- same.  In  other  words,  the  evaluation  will  consist  of 
comparing  the  sample  size  requirements  of  Weinstein's, 

Shanmugam  and  Balaban's,  and  the  Modified  Extrapolation 
techniques  for  a  given  variance  and  probability  of  error. 

For  clarity's  sake,  the  pure  Monte  Carlo  methods  will  also 
be  included  in  the  comparison. 

Parameters  and  constants  used  in  the  analysis  that  is  to 
follow  have  been  selected  on  the  basis  of  consistency  with 
those  selected  and  derived  by  the  authors  of  the  methods  to 

which  they  apply.  The  desired  threshold,  K,  was  chosen  to 

—  6 

yield  a  PM(e)  =  10  .  For  the  constant  variance,  it  was 

decided  that  the  selection  of  a  normalized  error  [see  equations 
(139)  and  (140)]  of  0.1  was  reasonable  as  well  as  computationally 
convenient.  The  bias  applied  to  the  uplink  noise  pdf  is 

B (x)  =  c/[fn(x)]a  (168) 

# 

where  c  is  a  constant  such  that 

c  =  [  ( 1-at )  /  (2 it )  a ]  ** .  (169) 

From  Shanmugam  and  Balaban,  values  for  x  and  c  have  been  chosen. 

The  analysis  itself  will  consist  of  two  parts.  To 
begin,  expressions  for  the  variances  of  each  of  the  estimators 
will  be  either  given  or  derived  to  an  extent  that  should 
provide  sufficient  justification.  The  second  part  will  be 
the  actual  specification  of  thresholds,  constants,  etc., 
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followed  by  the  calculation  of  the  sample  size  requirement 
of  each  estimator  for  the  prescribed  conditions. 


4.2.1 


Variance  Calculations 


A  BER  estimator  variance  analysis  must,  of  course, 
begin  with  the  derivation  of  the  variance  of  a  pure  Monte 
Carlo  estimator.  Assuming  that  the  estimator  is  to  be 
applied  to  a  binary  system,  as  has  been  done  throughout 
this  study,  then  the  threshold  test  on  each  received  bit  is 
an  obvious  Bernouli  trial.  To  illustrate  the  usefulness  of 
this  fact,  assume  that  the  data  to  be  transmitted  are  all 
0's.  Let  the  output  of  the  decision  logic  be  Y  =  y^  a  random 
variable  with  a  binomial  distribution.  The  expected  value 
and  the  variance  of  Y  are  then  known  to  be 


E[Y]  =  u v  =  np 


E[Y2]  -  E [Y] 2  =  a2  =  np(l-p) 


(170) 


(171) 


where  n  is  the  number  of  bits  transmitted  and  p  is  the  bit 
probability  of  error.  When  p  is  to  be  estimated  using  the 
error  counting  technique  with  N  samples  being  examined, 

C 

then  the  estimator  of  p,  p',  is 


p--i 


(172) 


I 

The  mean  of  this  estimator  is 


(173) 


and  the  variance  is  thus 


♦ 


(174) 


Since  the  probabilities  of  error  being  dealt  with  are 
PM(e)  <<  1,  the  variance  may  be  approximated  as 


The  variance  of  the  modified  Monte  Carlo  estimator 
of  Shanmugam  and  Balaban  is  quite  easily  approximated  by  the 
integral 


a .  % 
is 
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W  (x)  fx (x) dx 


(176) 


where  W(x)  =  bTxT  t*ie  wei?hting  function  for  the  output 
and  fx(x)  is  the  unbiased  pdf  of  the  input  noise.  Unfor¬ 
tunately,  as  personal  experience  and  conversations  with 
K.  S.  Shanmugam  bear  out,  this  expression  is  not  analytically 
tractable  in  general.  Shanmugam  and  Balaban,  however,  have 
generated  sufficient  data  to  make  a  plot  of  sample  size 
reduction  factor,  r,  versus  PM(e)  which  is  reproduced  here 
as  Figure  9 .  Since  the  determination  of  r  requires  a  fixed 
normalized  error,  e,  it  is  possible  to  determine  an  empirical 

expression  for  the  importance  sampling  variance. 

2 

The  procedure  taken  to  find  ar  ,  the  variance  of  the 
sample  size  reduction  factor,  is  essentially  just  that  needed 
to  determine  the  equation  of  the  line  on  the  plot  shown  in 
Figure  9.  Both  axes  are  seen  to  be  on  a  logarithmic  scale 
so  that  the  equation  of  the  line  is  log  y  =  m  log  x  +  b. 

The  solution  proceeds  as  follows.  It  is  known  that 
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y  =  r  = 


x  =  PM(e) 


Nc  = 
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Figure  9  -  Plot  of  Modified  Monte  Carlo  Simulation 
Results:  Sample  Size  Ratio  vs.  PM(e) 


Reproduced  from  K.  S.  Shanmugam  and  P.  Balaban,  "On  a  Modified 
Monte-Carlo  Simulation  Technique  for  Estimating 
Error  Rates  in  Digital  Communications  Systems," 
IEEE  Tr.  on  Commun. .  vol.  28,  Nov.  1980, 
pp.  1916-1921. 


At  this  point,  the  expression  for  the  importance  sampling 
variance  may  be  "derived”  via  the  step  shown  below. 


log  r  *  m  log  PM(e)  +  b 
r  =  (iob) [PM(e) ]m 


N 


c 


rNis 


Nis(1°b)  [PM(e)  ]I" 

Ve) 

Nc 

[pM(e) ]1_m 

lobNis • 


(177) 


From  Figure  7,  m  is  found  to  be  -0.6.  Selecting  the  point 
on  the  line  where  r  =  400  and  PM(e)  »  10~6,  b  is  found  to 
be  -1.0  and  the  variance  is  thereby  written 


10 


(PM(e)] 
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Variance  expressions  for  two  of  the  four  performance 
predictors  being  compared  in  this  analysis  have  been  determined. 
The  third  method  to  be  looked  at  is  Weinstein's  tail  probabili¬ 
ties  approach  [40],  To  begin  the  derivation,  recall  the 
three-threshold  case  and  procedures  described  in  the  previous 
chapter.  Note  that  yR  ■  ln(-ln [PM(e | K) ] )  where  K  is  the 


threshold  for  which  it  is  desired  to  estimate  the  BER. 


From  Weinstein's  paper  we  know  that  the  estimate  of  yR,  y',  i: 


m  aiY 


1J1 


+  a-y 


2J2  +  a3y3 


(179) 


where 


In  K  -  In  x 2 
In  x3  -  In  x^ 
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1 

a2  =  T 

,  In  K  -  In  x, 

»i  ■  T  *  In  x-— In  I"  “  2‘097  (180) 

The  variance  for  the  tail  extrapolation  method  is  easily 
seen  to  be 


var(yR)  =  |  j  a^  covfy^y..)  (181) 

Weinstein  notes  that  for  small  perturbations  in  P^(e|xj), 

Ayj  %  (dPM(eix.)]APM(e|xj) 

*  [P^elXj)  In  PM(e|xj)]"1  aP'telXj),  (182) 


If  the  covariance  of  the  BER  estimates  is  small,  it  is  seen 
that 


cov  [P^tejXj^)  ,P^(e  |xj)  1 


cov(yify.)  %  P  (e|x,)PM(e|x.)  In  Pm(e|xi)  In  P»(e|x.) 


(183) 
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PM^elxn^Nte 


PM(e lx±)  PM(®|Xj)  In  pM^e I xi)  ln  PM(e‘xj 


(184) 


where  n  =  max(i,j)  and  Nfce  is  the  tail  extrapolation  sample 
size.  With  the  assumption  that  the  (PM(e|xm)}  have  Gaussian 
distributions,  (184)  can  be  plugged  back  in  to  (181)  and  the 
variance  of  the  estimator  may  be  calculated.  Although  using 
(181)  promises  to  be  rather  difficult,  it  may  be  verified 
numerically  that  the  i  =  j  =  3  term  dominates  and  thus,  the 
variance  of  yR  may  be  written 


var(y  )  %  — — -  “7  N“PM(elX3) 

[PM(e|x3)  lnPM(e|x3)]  te 

=  f— -7 — | - r—j — =—7 — | - r)  var [P'  (e|x,)  ] 

^PM(e|x3)  ln  PM(e|x3) J  M  3 


(185) 


from  which  it  is  seen  that 


^PM(e|K)  ln  PM(e|K)  ^ 
\PM(e  I  X3}  ln  PM(e|X3)  / 


var[P^(e|x3) ] .  (186) 


The  final  variance  derivation  is  for  the  Modified 
Extrapolation  technique  that  has  been  proposed  in  this  chapter. 
Derivation,  however,  is  too  strong  of  a  word  for  what  needs 
to  be  done  to  determine  omg.  If  it  is  remembered  that  the 
proposal  is  to  use  the  Modified  Monte  Carlo  method  to  estimate 
the  BER  for  the  third  extrapolation  threshold  used  in  Weinstein's 
approach,  then  the  determination  becomes  the  substitution  of 
(178)  for  the  var (P^(e |x3) ]  in  (185).  The  expression  for  the 


variance  of  the  proposed  estimator  is  then 
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4.2.2  Performance  Comparisons 

A  Gaussian  distribution  is  assumed  from  which  PM(e)  is 

drawn  for  the  following  comparisons:  With  PM(e|K)  chosen 

to  be  1<T6,  K  -  4.78.  Since  it  has  already  been  stated  that 

105  samples  are  at  the  upper  limit  of  feasibility  for  pure 

Monte  Carlo  testing#  the  third  threshold#  x.,,  will  be 
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selected  to  yield  a  PM(ejx3)  »  10  :  x3  *  3.70.  Using  the 

requirement  that  x^#  x2#  and  x3  be  selected  to  be  evenly 
spaced  on  a  logarithmic  scale  along  with  what  is  known  from 
the  equations  for  the  a^s,  the  remaining  two  thresholds  are 
chosen  as  x^  =  3.02  and  x2  ■  3.35. 

A  normalized  error  of  0.1  has  already  been  specified 
so  that  the  sample  size  requirements  may  be  computed  as 
follows  from  one  of  the  two  definitions  of  e: 


2  -  1 
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For  the  error  counting  case# 


N 


e  P„(e  j  K) 


102 


Shanmugam  and  Balaban's  method  leads  to 


N.  =  10 

is  e5(PM(e|K)  4 

while  the  tail  extrapolation  approach  yields 

a32  /  In  PM(e|K)  \  2 

Nte  =  ^PM(e|x3)  m  PM(e|x3)j  PM(elx3)* 

Finally,  the  Modified  Extrapolation  technique  requires  a 
sample  size  of 

a,2  /  In  PM(e|K>  \  2 

"me  =  ^PM(e|x3)  In  PM<elx3)j  10[PM(«|x3>] 


Table  1  presents  the  results  of  these  calculations. 


TABLE  1 


COMPARISON  OF  SAMPLE  SIZE  REQUIREMENTS 


P  (s  K)  e.  Method  ♦  Samples 
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SUMMARY  AND  CONCLUSIONS 


* 
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In  the  past  four  chapters  a  fairly  thorough  study  of 
the  search  for  an  "ideal"  performance  prediction  tool  has 
been  made.  This  study  has  attempted  to  present  a  comprehensive 
view  of  the  need  for  an  efficient,  accurate,  in  short,  a 
near-ideal  predictor  for  use  in  the  analysis  and  design  of 
digital  communications  systems  both  in  general  and  in  specific 
application  to  satellites.  Chapters  1  and  2  laid  the  foundation 
for  an  understanding  of  the  basic  problem  that  confronts 
satellite  communications  systems  designers  and  the  need  for 
a  performance  prediction  tool  that  is  accurate,  flexible, 
easy  to  use,  and  computationally  efficient  (fast) .  In 
Chapter  3  a  fairly  comprehensive  view  of  the  current  state 
of  affairs  is  presented  along  with  a  review  of  past  work  and 
its  historical  development.  Finally,  in  Chapter  4,  the  best 
of  the  tools  that  have  already  been  proposed  are  identified 
and  a  combination  of  these  is  proposed  and  examined. 

It  was  discovered  that  there  are  two  distinct  categories 
that  the  extant  performance  prediction  tools  fall  into: 
the  analytical  approaches  and  the  simulation-based  approaches. 
While  the  former  were  recognized  as  fulfilling  the  speed 
requirements,  they  were  shown  to  be  inadequately  flexible 
and,  with  the  increasing  use  of  TDMA,  insufficiently  accurate 
vis'  the  effects  the  nonlinearly  operated  TWTAs  have  on  the 
noise-corrupted  uplink  signal.  On  the  other  hand,  the  simula¬ 
tion-based  approaches  were  shown  to  be  computationally  too 
slow  or  else  relatively  inaccurate  with  the  two  characteristics 
in  a  mutually  exclusive  trade-off.  Overall,  however,  the 
simulation-based  tools  were  found  to  be  superior  to  the 
analytical  ones  and  two  were  identified  as  fulfilling  the 

105 


greatest  number  of  the  requirements  presented  in  the  first 
two  chapters:  Weinstein's  tail  extrapolation  method  and 
Shanmugam  and  Balaban's  modified  Monte  Carlo  approach. 

In  an  attempt  to  create  a  new  and  better  performance 
predictor,  it  was  noticed  that  the  two  tools  identified  as 
having  the  most  potential  seemed  to  compliment  one  another. 
Consequently,  it  was  proposed  that  a  new  prediction  tool, 
called  the  Modified  Extrapolation  technique,  be  applied  using 
what  were  seen  as  the  best  of  the  characteristics  of  both 
Shanmugam  and  Balaban's  and  Weinstein's  tools.  The  proposed 
new  predictor  was  described  and  a  variance/sample  size 
reduction  analysis  was  performed  to  check  the  validity  of  the 
proposal.  It  was  discovered  that  the  Modified  Extrapolation 
method,  while  not  as  good  as  Shanmugam  and  Balaban's  approach, 
is  indeed  better  than  anything  else  that  has  previously 
existed  or  been  proposed.  However,  the  sample  size  reduction 
is  less  than  3dB  worse  than  the  modified  Monte  Carlo  technique 
for  the  normalized  errors  selected  and  examined  in  the  last 
chapter.  When  all  of  the  assumptions  that  have  been  made 
during  the  formulation  and  analysis  of  the  two  methods  are 
taken  into  account,  it  is  reasonable  to  state  that  the  difference 
in  required  sample  sizes  is  insignificant.  Due  to  the  increase 
in  complexity  presented  by  the  new  technique  over  its  close 
competitor,  one  must  conclude  that  the  modified  Monte  Carlo 
technique  of  Shanmugam  and  Balaban  is  superior  for  the 
normalized  error  and  error  probability  of  the  typical  satellite 
communications  system.  Remember,  however,  that  the  sample 
size  reduction  factors  of  both  Weinstein's  method  and 
Shanmugam  and  Balaban's  approach  increase  as  the  BER  decreases. 


It  can  be  seen,  therefore,  that  the  Modified  Extrapolation  will, 
for  some  BER  <  10“® ,  be  superior  to  the  importance  sampling 
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